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FOREWORD 


This  manual  is  a supplement  to  the  Engineering,  User's 
and  Systems  manuals  prepared  for  the  Antenna  Modeling  Program 
(AMP),  and  describes  the  operation,  theory  and  coding  of  the  changes 
made  to  AMP  for  more  accurate  treatment  of  multiple  wire  junctions 
and  reduction  of  the  time  for  interaction  calculations  on  large  structures. 

The  AMP  code  as  modified  (AMPJ)  has  been  delivered 
to  the  Naval  Ship  Engineering  Center  and  U.  S.  Army  Strategic  Communi- 
cations Command  and  was  developed  under  Office  of  Naval  Research 
Contract  N00014-71-C-0187. 
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1.0  INTRODUCTION 


Wire  antennas  and  their  supporting  structures  frequently  include 

junctions  of  several  wires  connected  together  at  a point.  Such  multiple 

junctions  are  especially  common  when  wire  grids  are  used  to  model 

solid  surfaces  as  is  often  done  in  mathematical  modeling  of  antennas. 

12  3 

The  antenna  modeling  computer  program  AMP  ' ' for  wire  antennas 

4 

and  the  extended  program  AMP2  for  modeling  wires  and  surfaces  both 
allow  for  modeling  multiple  wire  junctions,  and  have  demonstrated  good 
results  for  many  such  structures.  In  some  cases,  however,  when 
segments  of  greatly  different  lengths  have  been  joined  at  multiple  junctions 
the  method  used  in  these  programs  has  yielded  inaccurate  results.  Hence 
a more  careful  treatment  of  the  current  interpolation  at  multiple  junctions 
was  developed  and  is  included  as  an  option  in  program  AMPJ.  The 
details  of  this  method  are  given  in  Section  3.0  of  this  manual  and  the 
data  cards  to  request  the  optional  treatment  at  specific  junctions  are 
described  in  Section  2.  0.  Although  the  optional  technique  is  numerically 
more  stable  for  the  general  case  of  unequal  segment  lengths  at  a junction, 
it  was  not  included  in  the  extended  code  AMP2  for  the  following  reasons: 

1)  the  technique  used  in  AMP2  has  demonstrated  good  results  for  equal 
segment  lengths  and  has  a good  record;  2)  the  new  technique  as  imple- 
mented for  testing  requires  an  extra  unknown  at  each  multiple  wire 
junction  where  used,  and  there  are  many  junctions  in  a wire  grid;  3)  the 
new  technique  requires  more  development  and  testing  which  was  outside 
the  scope  of  this  contract. 

Program  AMPJ  also  includes  an  approximate  matrix  filling 
method  that  may  be  used  for  interactions  between  segments  separated 
by  more  than  a specified  distance  to  reduce  matrix  fill  time.  Hence 
the  only  feature  of  AMP2  that  is  not  available  in  AMPJ  is  modeling  of 

surfaces  via  the  magnetic  field  integral  equation.  AMPJ  includes  all 

% 

features  of  program  AMP  and  in  addition  the  optional  junction  inter- 
polation and  approximate  matrix  fill  method. 


Section  2.  0 of  this  manual  contains  instructions  for  use  of  the 

features  of  AMPJ  not  in  program  AMP,  and  supplements  the  AMP 

Users  Manual*.  Section  3.0  gives  the  equations  for  the  junction  inter 

polation  and  approximate  matrix  fill  methods,  and  supplements  the 

2 

AMP  Engineering  Manual1”.  Finally,  section  4.  0 details  the  coding  of 
the  routines  that  differ  from  those  in  program  AMP,  supplementing 

3 

the  AMP  Systems  Manual  . .Listings  of  the  changed  routines  are 
included. 


2.  0 PROGRAM  OPERATION 

The  basic  information  needed  to  use  program  AMPJ  is  contained 
in  the  AMP  Users  Manual*.  This  section  contains  supplementary 
instructions  and  information  for  using  the  optional  junction  interpolation 
and  approximate  matrix  fill  methods.  If  the  new  options  are  not 
required,  AMPJ  may  be  used  exactly  as  AMP.  The  one  exception  is 
that  AMPJ  uses  a time  saving  approximation  in  filling  the  interaction 
matrix  for  interaction  distances  greater  than  one  wavelength.  For 
results  identical  to  those  of  AMP  in  all  digits  printed,  this  approximation 
range  should  be  increased  to  greater  than  the  maximum  structure 
dimension  in  wavelengths  by  use  of  a KH  data  card. 

The  standard  interpolation  method  involves  the  extrapolation 
of  a segment  current  by  the  average  distance  of  the  centers  of  other 
segments  connected  to  the  junction,  requiring  that  the  extrapolated 
current  for  the  segment  equal  the  negative  of  the  sum  of  the  currents 
at  the  centers  of  the  other  segments.  Because  of  the  average,  this 
technique  can  lead  to  problems  when  the  segment  lengths  at  the  junction 
are  greatly  unequal.  The  optional  junction  interpolation  method  satisfies 
Kirchhoff's  Current  Law  directly  at  the  junction  and  also  forces  the 
derivatives  of  the  current  with  respect  to  distance  at  the  ends  of  each 
of  the  segments  at  the  junction  to  be  equal.  The  condition  on  current 
derivatives  is  based  on  continuity  of  potential  for  equal  segment  radii. 

Since  an  additional  unknown  is  required  in  the  matrix  equation 
for  each  junction  at  which  the  optional  junction  interpolation  method  is 
used,  this  method  should  be  used  only  when  required  for  accuracy.  The 
standard  interpolation  method  has  been  found  to  work  well  as  long  as 
the  connected  segments  at  multisegment  junctions  have  nearly  equal 
lengths.  Results  in  Section  3.0  show  an  error  of  about  15  percent  in 
the  input  admittance  of  an  antenna  when  the  segment  lengths  at  a junction 
near  the  source  differ  by  a factor  of  two.  The  accuracy  of  the  computed 
radiated  field  is  less  sensitive  to  the  segment  lengths  than  is  current. 

•In  general,  however,  segment  lengths  should  be  kept  within  a factor  of 
two  when  the  standard  interpolation  method  is  used. 
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With  the  optional  junction  interpolation  method  limited  testing 
has  shown  that  good  accuracy  may  be  obtained  with  segment  lengths  at 
a junction  differing  by  factors  of  8 to  10.  With  the  AMPJ  program  an 
antenna  model  may  conveniently  be  run  both  with  standard  and  optional 
interpolation  methods  to  test  the  need  for  the  optional  method.  When 
complex  structures  are  being  modeled,  it  is  quite  advantageous  to  be 
able  to  use  different  length  segments  at  junctions  with  confidence.  As 
a result  it  may  be  possible  to  reduce  the  number  of  segments  so  that 
the  matrix  size  and  the  running  time  are  decreased  in  spite  of  the 
additional  unknowns  for  multiple  junctions. 

The  program  execution  time  is  the  same  as  for  program  AMP 
except  for  differences  due  to  use  of  the  optional  junction  interpolation 
method  and  approximate  matrix  fill.  The  central  processor  time 
approximately  follows  the  formula 

T = Ak  (1  - 0.  7 Rw)  NZ/M  + B (N  + N.)3/M2 
+ CNeN2/M  + DkNjN 

where  N = number  of  segments  in  model 

M = number  of  degrees  of  symmetry 

N = number  of  different  excitations 

e 

= number  of  far  field  calculation  points 

N.  = number  of  junctions  at  which  optional  inter- 
polation method  is  used 

R = the  fraction  of  all  segment  pairs  for  which  the 

separation  is  greater  than  R0  where  R0  is  the 
limit  set  on  the  KH  card  for  change  over  to  the 
approximate  matrix  fill  method. 

k = 1 if  structure  is  in  free  space 

2 if  structure  is  over  ground 

and  A,  B,  C,  D are  proportionality  constants.  The  first  term  in  this 
equation  represents  the  time  for  matrix  filling;  the  second  term,  matrix 
factoring;  the  third  term,  solution  for  the  current  and  the  fourth  term. 
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calculation  of  the  far  fields.  Each  term  represents  only  the  dominant 
component  neglecting  terms  of  lower  order  in  N.  The  proportionality 
factors  depend  on  the  computer  system  on  which  the  program  is  run. 

To  give  an  idea  of  the  importance  of  the  terms,  the  factors  in  seconds 
for  a CDC  6600  computer  with  the  program  compiled  under  the  Run 
compiler  and  the  matrix  fitting  in  core  are  roughly 

A = 2 x 10"3,  B = 5 x 10"6,  C = 2 x 10"5,  D = 3 x 10"4 

If  an  antenna  is  analyzed  for  only  a single  excitation  and  the  far  field  is 
computed  at  a few  angles  the  execution  time  will  consist  almost  entirely 
of  the  time  to  fill  and  factor  the  interaction  matrix.  If  a number  of 
excitations  are  requested,  especially  for  out  of  core  solutions,  the 
time  to  solve  the  factored  matrix  equation  for  the  current  distribution 
can  become  significant,  and  if  a large  number  of  far  field  calculations 
are  requested  their  computation  time  must  be  considered. 

2.  1 NEW  INPUT  CARDS 

The  input  to  program  AM PJ  is  identical  to  that  for  program  AMP 
except  for  the  addition  of  data  cards  to  specify  junctions  at  which  the  optional 
interpolation  of  the  current  is  to  be  applied,  arid  the  separation  distance 
at  which  the  matrix  filling  changes  over  to  an  approximate  form.  Hence 
the  user  should  refer  to  the  AMP  Users  Manual*  for  the  basic  input 
data  structure.  If  only  the  basic  data  cards  are  used  the  interpolation  at 
all  junctions  will  be  the  same  as  that  in  program  AMP  and  approximate 
matrix  fill  will  be  used  for  segments  separated  by  more  than  one  wave- 
length. 

To  specify  matched  derivative  interpolation  at  one  of  more 
junctions  a one  must  be  punched  in  column  10  of  the  GE  card  at  the  end 
of  the  geometry  data.  This  card  may  be  followed  by  one  or  more  cards 
with  the  mnemonic  JX  to  specify  the  junctions  at  which  matched  derivative 
interpolation  is  desired  and  finally  a card  with  mnemnic  JE  to  indicate 
the  end  of  junction  specifications.  The  form  of  these  data  cards  is 
shown  below. 
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END  GEOMETRY  INPUT  (GE) 


CARD: 


5|  iOl  2C 


301  40| 501  60| 70|  80| 


C£  I II  I 12  Blink  Blink  I Blink  | Blink  Blink  I Blink  I Blink 


The  numbm  ilonf  tht  lop  rcfcr  u>  tht  Uw  tokann  in  Cl  eh  fitkl. 


The  function  of  this  card  is  the  same  as  in  program  AMP  except 
for  the  addition  of  the  integer  12.  If  12  is  equal  to  1 the  program  reads 
junction  specification  cards  following  the  GE  card.  If  12  is  blank  standard 
data  cards,  described  in  the  AMP  User's  Manual  are  expected  after  the 
GE  card. 


PURPOSE: 


to  specify  segment  junctions  (simple  or  multiple)  at 
which  the  optional  interpolation  is  to  be  used. 


CARD: 


J* 


» 


IV 


15 


IZ 


PARAMETERS: 

INTEGERS 
IX  - 


IY 


IZ 


r* 

30 

40 

50 

60 

70 

80 

Blank 

Blank 

Blank 

Blank 

Blank 

Blank 

Blank 

The  numbers  along  rhe  top  refer  to  the  last  column  in  each  field. 

II  I I I 


tag  number  of  one  segment  at  the  junction. 

Blank  or  zero  for  IX  implies  that  the  segment 

will  be  identified  by  the  absolute  segment 

number  in  the  next  location  (IY). 

equal  to  m,  specifies  the  m segment  of  the  set 

of  segments  with  tag  numbers  equal  to  IX.  If 

IX  is  zero  or  blank,  IY  is  the  absolute  segment 

number. 

specifies  the  end  of  the  segment  determined 
by  IX  and  IY.  IZ  equal  to  2 specifies  end  2 
of  the  segment  and  1 or  blank  specifies  end  1 
(reference  direction  for  current  is  toward 
end  2). 


NOTES: 

• Optional  interpolation  is  used  for  the  junction  at  the  specified 
end  of  the  specified  segment.  The  junction  may  be  either 
simple  or  multiple.  A JX  card  is  required  for  only  one 
segment  end  at  a junction  to  cause  use  of  optional  inter- 
polation for  all  segments  at  the  junction. 


i ' 


If  IX  and  IY  are  blank,  the  card  will  cause  optional 
interpolation  at  all  multiple  junctions.  Such  a JX  card 
must  occur  before  any  other  JX  cards  specifying  simple 
junctions.  For  optional  interpolation  at  all  multiple 
junctions  but  no  simple  junctions  the  JX  card  may  be 
omitted  with  only  a JE  card  used. 

Each  junction  at  which  optional  interpolation  is  used 
adds  an  additional  unknown  in  the  matrix  equation. 

All  segments  at  a junction  at  which  optional  interpolation 
is  used  should  have  equal  wire  radii. 

When  using  the  optional  interpolation  method  on  a 
structure  for  which  symmetry  is  used  in  the  solution, 
any  junctions  to  be  specified  in  the  first  symmetric 
section  must  be  specified  first  followed  by  the  same 
junctions  in  the  second  section  and  continuing  through  all 
sections  in  the  order  that  they  were  produced  by  reflection 
or  rotation.  In  addition,  the  junctions  must  be  specified 
in  the  same  order  in  each  section.  If  only  a JE  card  (no 
JX  card)  is  used  to  specify  all  multiple  segment  junctions 
these  rules  will  automatically  be  satisfied.  Use  of  a JX 
card  with  IX  and  IY  blank  to  specify  all  multiple  junctions 
followed  by  other  JX  cards  for  simple  junctions  is  not 
allowed  with  symmetry  since  the  junctions  will  not  be 
specified  in  the  proper  order. 
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PURPOSE:  to  mark  the  end  of  the  JX  cards  or  specify  all 

multiple  junctions. 


PARAMETERS:  None 

If  a 1 is  punched  in  column  10  of  the  GE  card  then  a JE  card  is 


required  to  return  the  program  to  reading  the  standard  input  cards. 

If  one  or  more  JX  cards  are  used  the  JE  card  is  placed  at  the  end  of 
the  JX  cards  and  its  only  function  is  to  mark  the  end  of  these  cards. 

If  no  JX  cards  are  used  the  JE  card  alone  will  cause  the  use  of  optional 
interpolation  at  all  multiple  junctions  (three  or  more  wires  joined). 


INTERACTION  APPROXIMATION  RANGE  (KH) 


PURPOSE:  to  set  the  minimum  segment  separation  distance  for 

U3e  of  a time  saving  approximation  in  filling  the  matrix. 


PARAMETERS: 

DECIMAL  NUMBERS 


RKH  - The  approximation  is  used  for  interaction 

between  two  segments  separated  by  more  than 
RKH  wavelengths. 

NOTES: 

• For  segments  separated  by  more  than  RKH  wavelengths 
the  interaction  field  is  computed  from  an  impulse  approx- 
imation to  the  segment  current.  The  field  of  a current 
element  located  at  the  segment  center  is  used.  For 
separations  less  than  RKH  a current  interpolation  function 
is  integrated  over  the  segment  length  as  in  the  basic 
AMP  program. 

• The  KH  card  can  be  placed  anywhere  in  the  data  cards 
following  the  geometry  and  junction  cards  (with  FR,  GN, 
EX,  etc.)  and  affects  all  calculations  requested  following 
its  occurrence.  The  value  of  RKH  may  be  changed  within 
a data  set  by  use  of  a new  KH  card. 

• If  no  KH  card  is  used  RKH  has  a default  value  of  1 wave- 
length. Hence  to  exactly  duplicate  a run  with  the  AMP 
program  a KH  card  should  be  used  with  RKH  greater 
than  the  maximum  structure  dimension. 


* 
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• The  minimum  value  of  RKH  which  can  be  used  to  obtain 
results  within  a few  percent  of  the  no  approximation 
case  seems  to  depend  to  some  extent  on  the  structure 
size,  type,  segmentation,  and  excitation.  Values  of 
.25  wavelengths  or  less  have  been  found  acceptable  for 
symmetrically  excited  structures  and  electrically  small 
wire  grids;  on  the  other  hand,  values  up  to  .5  wave- 
lengths have  been  required  for  very  asymmetrically  fed 
structures.  No  exact  guidelines  have  been  developed 
for  RKH;  therefore,  it  is  best  to  experiment  on  any  given 
problem  type  if  a minimum  value  is  desired.  RKH  should 
never  be  less  than  the  length  of  the  longest  segment, 
however. 

• The  matrix  fill  time  using  the  RUN  compiler  on  a CDC 

6600  computer  is  approximately  T = (2.  -1.4R  ) (10~^)  N^ 

1 w . 

seconds  where  R^  is  the  number  of  segment  pairs  for  which 

the  separation  is  greater  than  RKH,  divided  by  the  total 

2 

number  of  segment  pairs  (N  ).  Thus  the  fill  time  is 
decreased  by  about  70Rw  percent. 


t 
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2.2 


PROGRAM  OUTPUT 


The  program  output  is  essentially  unchanged  from  the  basic  deck. 
Segments  connected  to  junctions  at  which  optional  interpolation  has  been 
specified  are  indicated  in  the  block  of  segmentation  data  by  connection 
numbers  (1+  ot  I-)  less  than  -90,000.  Also,  the  value  of  RKH  is  printed 
following  the  printing  of  frequency  and  wavelength. 


! 


-12- 


3.  0 FORMULATION  OF  THE  1NTERPO.LATION  METHOD  AND 
APPROXIMATE  MATRIX  FILLING 

3.  1 INTERPOLATION  METHOD 

For  the  current  interpolation  method  used  in  program  AMP  the 
current  on  segment  j is  approximated  as 


I.  (s)  = A.  + B. 
3 3 J 


sin  k (s  - s.)  + C.  cos  k (s  - s.) 

i 3 3 


U) 


where  k = free  space  wave  number  (2 tt  / X ) 

Sj  = s at  the  center  of  segment  j. 

A.,  B.  and  C.  are  constants  to  be  determined  so  that  equation  (1)  yields 
3 3 3 

the  best  possible  approximation  to  the  true  current  on  the  segment.  Of 

the  3N  constants  to  be  determined  for  a structure  having  N segments, 

2N  are  eliminated  by  enforcing  conditions  on  the  local  behavior  of  the 

current.  These  conditions  are  used  to  eliminate  the  constant 

and  C.  in  terms  of  the  current  at  the  center  of  each  segment,  I.  = I.(s.). 

3 3 3 3 

The  N unknowns,  L,  are  then  computed  by  solution  of  the  matrix 

equation  derived  from  the  electric  field  integral  equation. 


In  the  basic  program  the  condition  used  to  eliminate  two  of  the 
three  unknowns  for  a segment  is  that  equation  (1),  when  extrapolated 
forward  over  the  distance  to  the  center  of  the  next  segment,  must  match 
the  current  of  the  center  of  that  segment,  and  when  extrapolated  back 
must  match  the  current  at  the  center  of  the  previous  segment.  At  a 
multiple  junction,  where  an  end  of  segment  j is  connected  to  two  or 
more  segments,  equation  (1)  is  extrapolated  a distance  equal  to  the 
average  of  the  distances  from  the  center  of  segment  j to  the  centers 
of  all  other  segments  connected  to  the  junction  and  required  to  equal 
the  algebraic  sum  of  the  currents  at  the  centers  of  the  other  segments, 
relative  to  the  reference  direction  of  segment  j.  These  conditions  are 
based  on  the  continuity  of  current  at  the  junction  (Kirchhoff's  Current 
Law).  Applying  these  conditions  to  equation  (1)  yields  the  three 
equations 


. ..  . 
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(2) 


A.  - B.  sin  k 57  + C.  cos  k 67  = K. 

3 3 3 3 3 3 

A.  + B.  sin  k fi.+  + C.  cos  k ft  = Kt 

3 3 3 3 3 3 


where 


A.  = half  the  length  of  segment  j 

rj*V=  ?*. 

V V°+  ? 4k 

k:=  ?(±v 

k%  E (±  i J 


The  summation  index  l takes  on  the  values  of  the  numbers  of  all  segments 

connected  to  the  first  or  - end  of  segment  j,  of  which  there  is  a total  of 

n , and  !<  takes  on  the  values  of  the  numbers  of  all  segments  connected 

to  the  second  or  + end  of  which  there  is  a total  of  n_^_ . The  plus  sign  in 

the  summation  of  currents  is  used  when  the  reference  directions  for 

segment  j and  segment  1 or  k are  parallel  and  the  minus  sign  when 

reference  directions  are  opposed.  Solving  equations  (2),  (3)  and  (4)  for 

A.,  B.  and  C.  yields 
3 3 3 1 

A.  = -7-  |*K7  sin  k it  - I.  sin  k (67  + fi+ ) + Kt  sin  k 6 . I (5) 

3 u L 3 33  33  3 3-1 

B.  = — |K7  (cos  k fit  -1)  + I.  (cos  k fi7  - cos  k fit) 

J A L j 3 3 3 3 

+ Kt  (1  - cos  k 6^  )J  (6) 

C.  = ~ |K7  sin  k fit  - I.  (sin  k fi7  + sin  k fit  ) + Kt  sin  k 6- 1 (7) 

3 - L 3 3 3 3 3 3 3j 

where  A-  sin  k fi7  + sin  k fit  - sin  k (fi^  + fit)  (8) 
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(9) 


f 

m 

■ 

■ 

I 

! 

m 


To  solve  for  the  L the  terms  are  regrouped  as 
I (s)  = K'  X (s)  + I Y (s)  + K+  Z (s) 

J J J J J J J 


Xj  (s)  = £sin  k 6 j + (cos  k 6*  -1)  sin  k (s  - s^)  - sin  k fit 
cos  k (s  - Sj)J 

(s)  = sin  k (6j  + 6^)  + (cos  k - cos  k 6j)  sin  k (s  - s^) 


+ (sin  k 6 j + sin  k 6 ) cos  k (s  - s 


;>] 


Z.  (s)  = — fsin  k 6 “ + (1  - cos  k 6’)  sin  k (s  - a.)  - sin  k 6" 

J A L 3 i 3 3 


(10) 


(11) 


(12) 


cos  k (s  - s 


•A 


The  electric  fields  for  filling  the  interaction  matrix  are  obtained  in  the  form 

of  equation  (9)  by  replacing  the  constant,  sin  k (s  - s.)  and  cos  k (s  - s.) 

3 3 

terms  in  equations  (10)  through  (12)  by  the  fields  at  the  observation  point 

due  to  these  current  distributions.  The  coefficient  of  each  1.  then  represents 

3 

a contribution  to  the  matrix  element  in  row  i and  column  j where  i is  the 
segment  at  which  the  field  is  evaluated. 

The  optional  interpolation  method  enforces  the  following  two 
conditions  at  a junction. 

1.  The  sum  of  currents  leaving  the  junction  is  zero. 

2.  The  derivatives  with  respect  to  distance  at  the  ends 
of  all  segments  at  the  junction  are  set  equal. 

These  conditions  are  applied  at  the  junction  rather  than  by  extrapolating 
to  the  segment  centers,  thus  eliminating  the  discontinuities  at  the  junctions. 
The  second  condition  is  based  on  the  continuity  of  potential  as  stated  in 
reference  5.  For  equal  wire  radii,  continuity  of  potential  implies  that 
the  charge  densities  on  the  wire  ends  at  a junction  are  equal,  which 
through  the  continuity  of  current  law. 


j-y  I (s)  = - jwq  (s) 


(13) 


implies  that  the  current  derivatives  are  equal. 


To  apply  these  two  conditions  at  a junction  an  additional  unknown 
is  introduced  representing  the  derivative  of  the  current  on  the  end  of 
each  wire  at  the  junction.  For  the  optional  interpolation  on  the  positive 
end  of  segment  j the  conditions  used  to  determine  A^,  and  in 
equation  (1)  are 

I.  (s.  - 6")  = K:  (14) 

3 3 3 3 

I.  (s.)  = I.  (15) 

J J J 

1.  (s.  + A.  ) = (16) 

3 3 3 3 

where  “ j = the  current  derivative  at  the  positive  end  of  segment  j 

The  solutions  for  A.,  B.  and  C.  are 

3 3 3 

A.  = [*K"  cos  k ^ . - I.  cos  k (A-  + 6 •)  + “t  s*n  ^ 6-  I (17) 

J D l J J J 3 J J J J 

B.  = 4=r  f(I.  - K~. ) sin  k A.  + (1  - cos  k 6 .)  1 (18) 

J D L J 3 3 3 3 J 

C.  = I.  - A . = f(I.  - KT)  cos  k A.  - “t  sin  k 6j"|  (19) 

3 3 3 D L 3 J J J Jj 

D - cos  k Aj  - cos  k (Aj  + 6j  ) (20) 

For  the  optional  interpolation  on  the  negative  end  of  segment  j the  conditions 
on  I(s)  are 

I . (s.  - A.  ) = oc~  (21) 

3 3 3 3 

I.  (s. ) = 1.  (22) 

3 3 3 

1.  (s.  + 6t)  = K+  (23) 

J J J J 

and  the  solutions  for  A.,  B.,  and  C.  are 

3 3 3 

A.  = 4=7  f«x  T sin  k 6 t + I.  cos  k (A  . + 6 t)  - K*  cos  k A-  I (24) 

3 D L J 33  333 

B^  = -jj  [«-  (cos  k 6j  - 1)  + (Ij  - K^)  sln  k 6 j J (25) 


For  the  optional  interpolation  on  both  ends  of  segment  j the  conditions  on 
I(s)  are 


and  the  solutions  for  A.,  B and  C.  are 


These  equations  give  the  constants  for  the  current  interpolation 
function  of  equation  (1)  once  the  current  values  L have  been  found.  They 
could  b'j  used  as  a starting  point  to  obtain  the  matrix  elements.  In  the 
program,  however,  the  matrix  elements  are  obtained  by  starting  with 
equations  in  the  form  of  equation  (9).  Taking  the  derivative  of  equation 
(9)  with  respect  to  s and  applying  the  appropriate  conditions  at  the  segment 
ends  leads  to  the  following  expressions  for  the  matrix  element  contri- 
butions for  the  row  corresponding  to  the  field  observation  point  on 
segment  i: 

1.  Optional  interpolation  on  + end,  standard  on  - end 


column  H I G.  X.  (s)  ds  - £ J G.  (s)  Z (s)  ds 
•'-A,  J-  J -Ai  1 J 

r A.  - A: 

column  j:  J * G.  (s)  Y-  (s)  ~ x~  J ^ ^ 

-A  1 J -A- 


column  oc 


: f * G.  (s)  Z.  (s)  ds  ‘ 

x 1 3 


2.  Optional  interpolation  on  - end,  standard  on  + end 

,A.  fA- 

column  k:  J * G.  (s)Z.  (s)  ds  - / G.  (s)  X.  (s)  ds 

-Aj  i J A 1 j 

A A- 

f j Y r ^ 

column  j:  G^  (s)  Yj  (s)  ds  - J ^ G^  (s)X^  (s)  ds 

A. 

column  «:  J G.  (s)  X.  (s)  ds 

X <.a  1 J 


3.  Optional  interpolation  on  both  - and  + ends 
A /"Aj 

column  j:  J * G.  (s)  Y.  (s)  ds  - J G.  (s)  X.  (s) 
-Aj  1 J . "Aj  1 J 

+ J G (s)  Z.  (s)  ds 

J A.  1 J 


column  ee  : 


X - Z 


A.  a 

/j  -A. 

A G (s)  X.  (s)  ds  - Z / J 
' '-A*  1 J •'-A 


Sj  (s)  ds) 


column  oc 


(s)  X.  (s)  ds  - 


For  each  case  G^  (s)  is  the  component  of  the  free  space  dyadic  Green's 
function  for  the  electric  field  tangent  to  segment  i at  the  center  of  that 
segment  due  to  a current  at  s on  segment  j.  X.  (s),  Y.  (s),  Z.  (s)  are 

J J J _ 

as  defined  in  equations  (10)  through  (12)  with  the  exception  that  6 • = A;  • 

J J 

when  optional  interpolation  is  used  on  the  - end  of  segment  j and 

fit  = A • when  the  optional  interpolation  is  used  on  the  + end,  and  s^  = 0. 

Also 

X = — - [cos  k (&.  + 6 ) - cos  k A . ] (43 

A l J j J J 

Y = — - cos  k (Aj  + 6 j ) J (44 


=-  r 

A l 


cos  k A- 


where  A is  defined  in  equation  (8)  and  6.  = 6.  if  optional  interpolation  is 

J J 

on  the  - end  only,  6j  = 6 j if  optional  interpolation  is  on  the  + end  only, 
and  6.  = A-  if  optional  interpolation  is  on  both  ends.  The  column  indicies 
l and  k take  the  values  of  the  numbers  of  all  segments  connected  to  the  - and 
+ ends,  respectively,  when  the  standard  interpolation  end  of  a segment  is 
a multiple  junction. 

The  columns  designated  by  « represent  the  unknown  current 
derivative  common  to  all  segment  ends  at  the  junction.  An  additional 
equation  for  this  unknown  is  obtained  from  the  derivative  of  equation  (9) 
evaluated  at  the  segment  end  at  which  the  current  derivative  is  «.  K ^ 
or  K".  is  replaced  by  the  current  at  the  segment  end.  By  forming  the 
sum  of  these  equations  for  each  segment  at  the  junction  the  currents  at 
the  segment  ends  are  eliminated  by  the  condition  that  their  sum  be  zero, 
leading  to  the  equation 
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1 


I 

j = 1 


© I Z/X  + I.  Y/X  + oc/x 

j J 


© - I,  Z/X  - I.  Y/X  + cc/X 

K . j 

J J 


© - L - ccX/(Z' 


z/(z‘ 


© I.  - « X/(Z* 
J 


X2)  - cc+  z/(z2  - x2) 


where  j = 1, M represents  the  numbers  of  all  segments  connected  to 

the  junction  and  the  expressions  apply  as  follows: 

® for  junction  = on  the  + end  of  segment  j,  standard 
interpolation  on  the  - end. 

© for  junction  cc  on  the  - end  of  segment  j,  standard 
interpolation  on  the  + end. 

® for  junction  « on  the  + end  of  segment  j,  and  junction 
K’  on  the  - end. 

© for  junction  = on  the  - end  of  segment  j,  and  junction 
«+  on  the  + end. 

The  above  equations  are  used  in  the  program  AMPJ.  It  is 
possible  to  eliminate  the  current  derivative  unknowns  («)  rather  than 
include  them  in  the  set  of  equations  solved  numerically  but  this  is  diffi- 
cult if  the  new  interpolation  method  is  used  on  both  ends  of  a segment. 
Also,  the  code  could  be  generalized  for  unequal  wire  radii  at  junctions 
but  this  has  not  been  done. 

The  stability  of  the  two  interpolation  methods  for  modeling 
multiple  junctions  is  shown  in  Tables  1 through  3 for  an  antenna  composed 
of  a linear  element,  fed  at  the  center,  with  vee  shaped  loads  on  each  end. 
The  input  impedance  computed  with  standard  interpolation  at  all  segment 
junctions  is  shown  in  Table  1 for  varying  segment  lengths  on  the  entire 
center  section.  For  this  model  there  are  4 segments  on  each  of  the 
arms  and  n segments  on  half  of  the  center  element.  Though  no  attempt 
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NUMERICAL  TEST  ANTENNA 


2 METERS 


>4 

METER 


WHEN  n = 1,  JUNCTION  SEGMENTS  ARE  APPROXIMATELY  EQUAL  IN  LENGTH 
2 SEGMENT  SOURCE  IS  USED 
TABLE  OF  INPUT  IMPEDANCE/2. 

TREQ.  (MHz) 


270 

280 

290 

1 

33.5  + J79.5 

36.7  + j95.4 

40.0  + jll  1. 

2 

35.8  + J93. 

39.5  + jl  1 0. 

43.4 + j1 26.8 

4 

33.2  + jll  9.6 

36.8  + jl  35. 

40.6+150.4 

6 

28.2 + j1 52.5 

31.1  + jl  65 . 

34.2  + j1 77.5 

8 

21.7  + jl90.2 

23.8  + j1 99.5 

26.0  + j208.6 

10 

13.9  + j233. 

15.1  + 1*237.9 

16.3  + j243.2 

TABLE  1 

INPUT  IMPEDANCE  OF  TEST  ANTENNA  USING  STANDARD  JUNCTION 
INTERPOLATION  SCHEME  FOR  VARIOUS  SEGMENT  LENGTHS  AT  JUNCTION 
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0425-13269 


1 METER 


has  been  made  to  keep  the  source  width  constant,  the  input  impedance  of 
the  structure  should  not  vary  as  rapidly  as  indicated  by  these  results. 


For  Table  2,  the  segment  on  the  center  wire  connected  to  the 
junction  is  set  equal  to  the  other  segments  at  the  junction  while  the 
remaining  center  wire  segments  are  progressively  decreased  in  length 
as  before.  In  comparison  with  the  results  in  Table  1,  the  results  in 
Table  2 are  quite  stable  to  a large  segment  length  factor.  This  serves 
to  validate  the  fact  that  the  junction  is  indeed  causing  the  problem,  and 
in  addition  it  shows  the  usefulness  of  the  standard  junction  interpolation 
technique  for  the  case  of  equal  segment  lengths  at  a multiple  junction 
and  unequal  segments  elsewhere. 

Table  3 contains  the  results  obtained  using  the  new  interpolation 
technique  at  the  two  multiple  junctions  and  the  standard  interpolation 
at  all  junctions  of  two  segments.  The  segments  on  the  center  section 
are  progressively  decreased  in  length  including  the  segment  connected 
to  the  junction.  These  results  show  much  greater  stability  than  the  results 
given  in  Table  1 where  the  segment  lengths  at  the  junction  are  varied  in 
a similar  manner.  On  the  other  hand,  the  results  for  the  new  technique 

i| 

are  somewhat  less  stable  than  for  the  case  of  equal  segments  at  the  junction 
(Table  2),  but  when  complicated  structures  are  being  modeled,  it  is 
quite  advantageous  to  be  able  to  use  different  length  segments  at  junctions 
with  confidence. 

3.  2 APPROXIMATE  MATRIX  ELEMENTS 

When  wire  segments  in  a structure  are  distant  from  an  obser- 
vation point  with  respect  to  wavelength,  simple  expressions  can  be  used 
to  obtain  accurate  values  for  the  fields.  This  fact  can  be  used  to  substan- 
tially reduce  the  time  required  in  calculating  the  corresponding  interaction 
matrix  elements.  The  following  expressions  are  used  in  the  AMPJ 
code  when  segment-observation  point  separation  permits: 
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TEST  ANTENNA 


SEGMENTS  AT  THE  MULTIPLE  WIRE  JUNCTION 
ARE  APPROXIMATELY  EQUAL  IN  LENGTH 


2 SEGMENT  SOURCE 


TABLE  OF  INPUT  IMPEDANCE/2. 
FREQ.  (MHz) *- 


270 

280 

290 

1 

34.4  + jS0.3 

37.9  + j96.7 

41.7  + jl  12.9 

2 

35.1  + j80.8 

38.8  + J97.5 

42.9 + j1 14.1 

4 

35.4  + 181. 

39.4  + j98. 

43.7+  j115. 

6 

35.6  + j’82.2 

39.5  + j98.6 

43.9  + j1 15.7 

8 

35.6+  j83.8 

39.6  + j98.6 

43.9  + ]1 12.4 

10 

35.9  + j99. 

39.6  + j97.8 

44.1  +j121.9 

. ■ » » 

TABLE  2 ‘ 7 

INPUT  IMPEDANCE  OF  TEST  ANTENNA  USING  STANDARD  JUNCTION 
INTERPOLATION  SCHEME  FOR  AN  EQUAL  SEGMENT  LENGTH  JUNCTION 


0425-13267 


TEST  ANTENNA 


s 


WHEN  n = 2,  JUNCTION  SEGMENTS  LENGTHS  ARE  APPROXIMATELY  EQUAL. 

2 SEGMENT  SOURCE 


CURRENTS  AND  DERIVATIVES  MATCHED  AT  JUNCTION 


TABLE  OF  INPUT  IMPEDANCE  /2. 


n 

1 

f = 270  MHi 

2 

34.5  + j78.9 

4 

36.1  + j77.6 

6 

37.0  + J76.9 

8 

37.7  + J76.1 

10 

38.1  + j75.4 

TABLE  3 

INPUT  IMPEDANCE  OF  TEST  ANTENNA  USING  THE  NEW  JUNCTION 
TECHNIQUE  FOR  UNEQUAL  JUNCTION  SEGMENT  LENGTHS 
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0425-13268 


Hgjfcgi 


1 


cos  0 


(47) 


Er(r0) 

E0fc>) 


1 1 

-ikr0 

( JQ_  + 

2 TT 

0 

L2 

Ii 

P*ikrt>  / 

' Ioju 

+ 

4 TT 

e \ 

iOJ  e r0 


ia)er0 


) sin  0 


(48) 


These  are  the  fields  of  an  incremental  dipole  of  moment  II  located  at  the 
origin  of  a standard  spherical  coordinate  system  and  oriented  in  the 
z direction^.  At  sufficient  distances  equations  (47)  and  (48)  are  used 
for  the  field  of  a segment  where  i is  set  equal  to  the  segment  length  and 
I is  set  equal  to  the  center  point  current.  Thus,  these  expressions  are 
the  same  as  would  be  obtained  using  a pulse  function  current  expansion 
and  one  step  integration. 

This  approximation  has  been  found  to  yield  good  results  for 
separation  distances  as  small  as  .25  to  .2  wavelengths.  Table  4 shows 
the  accuracy  obtained  for  a particular  structure,  a 2\  dipole,  for  various 
segmentations  and  for  various  separation  distances  for  which  the  expres- 
sions in  equations  (47)  and  (48)  were  used.  The  KH  parameter  in  the 
table  specifies  the  distance  at  which  change  over  to  the  approximate  field 
expressions  occurs.  The  column  on  the  left  hand  side  of  the  table  shows 
the  number  of  segments  away  from  the  field  segment  which  are  integrated 
over.  For  this  example  it  can  be  seen  that  the  impedance  accuracy 
remains  within  a few  percent  for  a KH  down  to  .21  wavelengths.  It  should 
be  pointed  out,  however,  that  due  to  the  quantized  nature  of  the  problem 
a KH  parameter  slightly  less  than  . 2 wavelengths  will  cause  an  abrupt 
change  to  integration  over  one  fewer  segments.  For  the  case  of  . Z\ 
segment  lengths,  this  means  integration  for  the  self  term  only  and  the 
results  are  poor.  This  problem  can  be  avoided  by  keeping  the  KH  para- 
meter larger  than  the  longest  segment.  It  should  also  be  pointed  out 
that  the  minimum  value  for  KH  seems  to  depend  to  some  extent  on  the 
structure  size,  type,  segmentation,  and  excitation.  Values  of  KH  up  to 
• 5X.  have  been  necessary  to  obtain  only  a few  percent  error  for  some 
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STRUCTURE  SEGMENT  LENGTHS 

NUMBER  OF 
SEGMENTS 
INCLUDED 

0.2 

0.1 

0.05 

0 

KH  PARAMETER 
.01 

.01 

.01 

% ERROR  REAL.  IMAG. 

47.2,  53. 

77.3,  135. 

97.4,  170.4 

1 

.21 

.11 

.06 

2.2,  .62 

12.4,  12.4 

21.4,  190 

.41 

.21 

.11 

2 

.068,  .015 

1.3,  2.4 

.12,31.5 

.61 

.31 

.16  j 

3 

.35,  .01 5 

.09,  .09 

.23,  9.5 

.81 

.41 

.21 

4 

•19..33 

.028,  .30 

.13,  3. 

6 

1.21 

.06,  .003 

.61 

.035,  .12 

.31 

.01,. 23 

8 

1.65 

.81 

.41 

.02,  .022 

.09,  .19 

.037,  .041 

Table  4 

PER  CENT  ERROR  OF  THE  INPUT  IMPEDANCE  OF  A 2X  DIPOLE 
USING  PARTIAL  INTEGRATION  AS  COMPARED  TO  COMPLETE 
INTEGRATION 


structures  with  very  asymmetric  feeds.  No  exact  guidelines  have  been 
established;  therefore,  it  is  probably  best  to  experiment  with  any  given 
class  of  problems  if  a minimum  value  of  KH  is  sought.  The  default 
value  for  the  KH  parameter  in  the  AMJPJ  code  is  one  wavelength. 


4.  0 DESCRIPTION  OF  COMPUTER  CODE  CHANGES 

The  following  subroutines  have  been  changed  to  implement  the 
optional  junction  interpolation  and  approximate  matrix  filling: 

CABC 

CMSET 

DATAGN 

FACTR 

INTG 

JMELS 

LFACTR 

MAIN 

SOLVES 

TRIO 

In  addition  the  variable  JMAX  has  been  added  to  common  block/DATA/ 
throughout  the  program.  JMAX  is  the  total  number  of  junctions  at  which 
the  new  interpolation  method  is  used. 

The  matrix  filled  by  subroutine  CMSET  consists  of,  first,  segment 
field  equations  in  the  order  of  segment  numbers  and  then  equations  for  the 
current  derivatives  from  equation  (46).  For  a structure  with  symmetry 
the  field  equations  for  the  first  section  are  followed  by  the  current 
derivative  equations  for  that  section  and  the  equations  continue  in  that 
order  through  all  symmetric  sections.  Thus,  for  a structure  with  two 
symmetric  sections  the  matrix  equation  has  the  form 


where  Ij  and  represent  colui..  vectors  of  the  unknown  currents 
and  current  derivatives,  respectively,  on  the  first  section.  Only  the 
upper  half  of  the  matrix,  representing  equations  for  the  first  section, 
is  stored  and  is  stored  in  transposed  form  as  in  program  AMP.  Sub- 
routine ETMNS  which  fills  the  right  hand  side  vector  has  not  been 
modified  and  hence  fills  the  applied  field  values  in  consecutive  locations. 
Subroutine  SOLVES,  however,  has  been  modified  to  insert  zeros  for 
the  current  derivative  equations  and  reposition  the  applied  field  values 
for  structures  with  symmetry  before  solving  the  matrix  equation. 

Some  of  the  common  block  lengths  have  been  changed  from 
those  in  program  AMP.  The  maximum  number  of  segments  is  800  in 
AMPJ  although  this  is  also  the  upper  limit  for  the  sum  of  the  number 
of  segments  and  the  number  of  junctions  at  which  the  optional  inter- 
polation is  used.  Also,  the  area  in  core  for  storage  of  the  interaction 
matrix  has  been  increased  to  10000  complex  numbers.  This  allows 
structures  with  up  to  100  unknowns  (segments  plus  junctions  with  new 
interpolation  method)  to  run  in  core.  In  this  form  the  program  requires 
approximately  240000g  words  of  storage  to  load  on  a CDC  6000  series 
computer  when  compiled  with  the  FTN  compiler. 

The  following  are  brief  descriptions  of  the  changes  to  the 
modified  subroutines.  Lists  of  these  routines  are  included  at  the 
end  of  this  section.  Since  the  routines  have  not  been  sequence  numbered 
the  changes  can  be  located  by  the  gaps  in  the  old  sequence  numbers. 
References  to  statement  labels  refers  to  the  labels  in  the  left  hand 
column  of  the  list. 


CABC 

This  routine  computes  the  constants  A^,  IL  and  CL  for  equation 
(1)  for  either  the  old  or  new  interpolation  method.  The  statements 
down  to  label  15+1  set  CLO,  CLL  and  CLY  as  follows: 


CLO  = K7  or  « ” 
J 

CLL  = I. 

J 

CLY  = K+.  or  « + 
J 
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A connection  number  less  than  -90000  for  a segment  end  indicates  the 
new  interpolation  method  is  used  at  that  end  requiring  or 
The  statements' from  label  16  through  21  compute  AX  = A^,  BX  = and 
CX  = Cj.  Statements  labeled  16  + 2 through  16  + 5 evaluate  equations 
(3)  and  (6);  17  + 1 through  17  + 4 evaluate  equations  (24)  and  (25); 

18  through  18  + 3 evaluate  equations  (17)  and  (18);  and  19  through  19  + 2 
evaluate  equations  (31)  and  (32).  For  each  case  statement  20  evaluates 
C.  = I.  - Aj.  Finally  from  statement  20  + 1 through  21  the  real  and 
imaginary  parts  of  the  constants  are  stored  in  arrays. 


CMSET 

Sequence  number 
references 

CM19.1+2  : 

CM19.  1+2  : 

CM122.  1 : 

CM127+1  : 


JSEQ(J)  = matrix  row  for  segment  J 

JJEQ(J)  = matrix  row  for  current  derivative 
at  junction  J. 

Branch  to  section  for  approximate  matrix  fill. 


/«. 

* G.  (s)  X.  (s)  ds 

j 

• • 4 f ^ 

CE2  = J Gi  (s)  Y.  (s)  ds 

CE3  = / ^ G.  (s)Z.  (s)  d3 
-A.  1 3 


CM127  + 11  to 
CM  127  + 13  ; 

label  24  + 1 to 
label  24  + 3 : 

label  26  to 
label  26  + 3 : 


expressions  (37),  (38)  and  (39)  are  evaluated 

expressions  (34),  (35)  and  (36)  are  evaluated 

equations  (40),  (41)  and  (42)  are  evaluated.  In 
the  above  three  cases  the  contributions  to  « 
columns  are  entered  into  the  matrix.  Other  contri- 
butions are  entered  in  the  following  code. 


1 


-3L 


I 
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label  27  to 
label  28 


label  29 

label  29  + 1 
to  label  30 


CM  144.  2 to 
CM145 

label  18  + 1 
to  label  66 


DATAGN 


fill  matrix  elements  for  - end  of  segment  (currents 
in  K") 

J 

fill  matrix  element  for  I. 

J 

fill  matrix  elements  for  + end  of  segment  (currents 

in  K ) 

J 


Approximate  matrix  fill  section.  Equations  (47) 
and  (48)  are  evaluated. 

Equation  (46)  is  evaluated  for  each  segment. 
JCAS  corresponds  to  the  cases  1,  2,  3 and  4 in 
the  equation. 


The  coding  from  statement  label  9+2  through  the  end  of  this 
routine  sets  the  connection  numbers  for  segment  ends  at  which  the  new 
interpolation  will  be  used.  The  JX  or  JE  data  cards  are  read  at  statement 
111.  Statements  21  + 1 through  207  set  the  connection  number  for  a 
segment  end  specified  by  a JX  card  and  for  all  other  segments  connected 
to  that  segment  end.  Statements  208  through  210  search  for  all  multiple 
junctions  and  reset  the  connection  numbers  for  the  new  interpolation. 

The  latter  section  is  entered  when  a blank  JX  card  occurs  first  or  when 
there  is  no  JX  card  but  only  a JE  card. 

The  variable  JMAX  is  used  to  count  the  number  of  junctions  at 
which  the  new  interpolation  is  specified.  The  connection  numbers  for 
all  segments  connected  to  junction  number  JMAX  are  set  to  -(90000  + 
JMAX).  On  exit  from  the  routine,  JMAX  is  left  as  the  final  number  of 
junctions  with  the  new  interpolation  and  passed  through  common/DATA/  to 
other  routines. 

FACTR 


Minor  modifications  have  been  made  to  FACTR  following  sequence 
number  FA31  and  at  FA59  to  use  temporary  variables  to  avoid  unnecessary 
evaluations  of  subscript  references. 
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INTG 


Statements  added  following  IG42  evaluate  X,  Y and  Z of  equations 
(43),  (44)  and  (45)  and  store  them  in  XM,  YM,  and  ZM  respectively. 

JMELS 

The  calculations  of  the  matrix  row  indicies  JPJ  and  JMJ  have 
been  changed.  While  JP(J)  and  JM(J)  represent  segment  numbers  JPJ 
and  JMJ  are  the  locations  in  the  matrix  corresponding  to  these  segments, 
taking  into  account  the  additional  matrix  rows  for  current  derivative 
unknowns. 

LFACTR 

Minor  modifications  have  been  made  to  LFACTR  at  LF58  and 
L.F92  to  use  temporary  variables  to  avoid  unnecessary  evaluations  of 
subscript  references. 

MAIN 

Sequence  number 
references 

MA64  to  MA65  : 


MAI  39  + 1 

MA188  to 
MA189 

SOLVES 

Statements  between  sequence  number  SS11  and  SSI  3 insert  zeros 
in  the  B vector  for  the  right  hand  side  of  the  matrix  equation  in  locations 
corresponding  to  the  current  derivative  equations.  For  structures  having 
symmetry  the  applied  field  values  are  relocated,  using  Y as  scratch 
storage,  to  make  room  for  the  zeros  within  the  vector. 


JPMAX  = number  of  junctions  with  new  interpolation 
in  one  symmetric  section 

NEQ  = total  number  of  unknowns 

NPEQ  = number  of  unknowns  for  one  symmetric 
section. 

RKH  = default  value  for  separation  distance  at 

which  matrix  filling  changes  over  to  approx- 
imate form. 

Set  new  value  for  RKH 
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Statements  between  sequence  numbers  SS67  and  SS68  rearrange  the 
solution  vector  for  the  case  of  a symmetric  structure  which  uses  the  new 
interpolation  on  some  junctions.  The  solution  vector  at  SS67  consists  of 
the  currents  for  the  first  symmetric  section  followed  by  the  current 
derivatives  for  junctions  on  the  first  section,  then  currents  for  the  next 
section,  and  continues  in  this  order  through  all  sections.  These  statements 
put  the  currents  in  consecutive  locations  so  that  L is  in  location  j,  with 
the  current  derivatives  in  consecutive  locations  following  the  last  current 
value. 

TRIO 

The  statements  following  TR12  andTRlS  have  been  added  so  that 
QIL  or  DIK  are  set  to  Aj  for  new  interpola: .or.  on  the  - or  + segment  end 
respectively. 

Lusts  of  the  changed  routines  follow. 
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C 
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15 
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If <NEQ.NE.N9C JIGO  TO  * 

If  UST.CQ.lt  GO  TO  7 
NC  x *N»f 3 

IC*SE»1 
go  ro  u 
nCCL«?*nP®L* 

ICASE*3 
GO  TO  !1 

If  Usr.EO.lt  GO  TO  9 
NC0L*NP*a 

sf ocs*s»eo 
IC4SC*? 

GO  TO  11 

C*Lt.  f^LOCH  lNaLSY*.HOSTM*NLSrM,l*ES9*.NP£G.S«>F0*  »Nl  » 
nC0L*2*hP9L* 

If  TlNT.EO.lt  GO  TO  10 
NC DLS*N°EO 
I CASE** 

GO  TO  11 
*.C0tS*?*NPSY»« 

!Ca$C*s 
CCsT  I N»*E 
N?L0*  x*N9L0KS 
NP9LK<-»NPRlK 
NL ASTA*NLAST 
NWOSXasRON 
NCOLX*sCOL 
PRINT  135 

FIlE  PREPARATION  FOR  OOT-OF-COWF  MATRIX  SOLUTION.  FILES  ■'EmOUnO 
ANO  ENOflLE  mRITIEs.  ANO  TESTING  FOR  RESTART 

If  USTART.E0.9t  GO  To  12 
CALL  uncat 

00  TO  14 

If  UCASE.LT.3I  00  TO  14 

00  13  1*1*7 
NUN I T * I T AP  ( I > 

®E« INO  NUNIT 
CSO  Fit  E NUNfT 
»£#IN0  NUNIT 
CONTINUE 
CONTINUE 

OEFAULT  VALUES  F 0®  INOUT  PARAMETERS  AND  FLAGS 

1 GO*  l 

Fmm2S*100. 
f— 2*300. 
nf  ro* 1 

PKM«1 . 

NLOAO*fl 

■SVMPs] 

IirrPaft 
NE  T«0 
NRAQL*G 
NE APa-1 
IPTFLGa-2 
IF APs-1 

7PArf*CMPLX(I..O.) 

ipr»r*^ 

l3£D*9 

PAIN  INPUT  SECTION  - STANDARD  »EAO  STATEMENT.-  Ju-PS  To  APPRO- 
PRIATE SECTION  FOR  SPECIFIC  PARAMETER  SET  u'»  . 

READ  115*  AIn.IT-PI UT-P2. ITmpi.ITMP4.Tmp1 .TMP2•T-P3*TNPN,TMP6•T- 
1P6 

MPCNT»"RCNT*l 

PRINT  I 37.  MPCNT.AlN.ITMPl .ITMR2.TTMP3.IT-R4.T-PI .T-P2.TmP3.T-pm. 
1 T-PS.T-P6 

If  (AlN.EQ.ATSTUn  GO  TO  17 

If  tAIs.EO.ATST ( 3| > Go  TO  1® 

If  tAIs.EO.ATSTt-M  GO  TO  ?2 

If  t AIs.EO.ATSTtSt I GO  TO  25 

If  (A/n.E0.ATST|9>>  GO  TO  29 

If  lAlN.EO.ATSTU-t!  GO  TO  29 

If  CAIn.EO.ATSTUSII  SO  TO  32 

If  !A!n.E0.ATST(7»#  GO  TO  39 

If  t AJn.EO.ATSTOM  Go  TO  33 

If  t AIs.EO.ATSTlRl » GO  TO  35 

If  tAlN.EO.ATSTtlOtt  GO  TO  37 

If (AlN.E3.ATST (17) >G0  TO  20l 
If  t AlN.EO.ATSTtl?! ) GO  TO  1 

If  (AIN.NE.ATSTC13U  GO  TO  15 

if  ut-pi.ne.oi  call  catlog 

STOP 

PRINT  13® 
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If  (NfPO.EO.O*  NFR0*1 

f — —2*T— P I 
OELfROaTMP? 

If  UPfO.EO.lt  7PnOR“^0. 

I G0*1 

iflom*? 

GO  TO  15 


Ml 

51 

|*  X 

9? 

m A 

93 

— A 

R« 

— A 

ss 

— A 

96 

— A 

97 

Ml 

59 

-A 

59 

-A 

1 GO 

MA 

131 

-A 

102 

M| 

103 

Ml 

1 94 

m a 

195 
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200 

I 


36 


IFlO#** 

|F  < | GO#GT .3)  100*1 
If  I c.0  TO  l* 

NE  T »Nr  T • 1 

IF  <NCT .LE.SE TMIJ  GO  TO  31 
PB  1ST  1*? 

STOP 

NT  *P(nf» » • £ 

T F (AIn.EO.AT'T  1*01  NTfP<NFl>*» 

IS*G1  l\£T>  *tSFGN0T  I ? hp  1 > I t*"V> 

I 5F  02  T T * I SF  GSO  I IT  »»  ) « I I *°4  » 

TN»(Vf  T»*T“P1 
T I I ! I S*  T ) * T h»2 
V 12»  IN*  T ) *ThB  1 
»|?I(NU|  »T*p<. 

V22P(N*T  > *T*DS 
V2FI  TN*  T I *T**PA 

IF  CNTtO|sFT).E0.l.0».TMPl.GT.'>.)  f.Q  |0  IS 
NT  *P  T NF  T ) * 3 
V!  1PINFT)  *-T**Pi 
GO  TO  IS 

PPINT  FONTPOL 

IPTFLC.*!T“P| 

I P T A(«*  I T M®  P 
1 PT  A&f  * I T**®! 

JP  T A6T * I T **p% 

IF  < ITMO<..£0.ni  IPf4GT  = IPTP0F 
GO  TO  IS 

NE*P  field  CALCULATION  PA9AHEI EOS 

if  ».not.tiflo#.Eo.«.ano.nfmo.nf.i» » GO  10  1* 

PPINT  143 

PPINT  144.  ISECNIJ) 

NEAP*  I THPl 
NP»  * I T*P2 
NOT* I Tmo3 
NBZ*IT*P<* 
xnp*t«40( 

VNP*T*40? 

7N3*T«p  3 

DinQzTvpl 

QYSOsTMPS 

0Zn9*Tmp6 
IFlO»=« 

IF  (NFO0.NE.lT  GO  TO  IS 
GO  TO  <42.47.54,72,73)  . I GO 

G»OUNO  pepoesentation 

IF  T.NOT. < IFLOW.EO.P.ANO.NFPO.Nr. 1) > GO  TO  3 6 
POINT  !<**•.  ISECNI2) 

EPS02*T**P1 
SIG2-T-P2 
CL  T *TMP  3 
CmJ*Tmp^ 

IFlO**9 
GO  TO  IS 

C STANOAQO  OSSEPVATION  ANGLE  PAPa^ETEOS 

C 

37  tfapbItmpi 

nth*i  r*«P2 

NPh*I Tm03 

IF  (NTh.EO.OT  NTH*  1 

IF  (nPn.EO.OT  NOH*I 

IP0*ITmO4/10 

I A*P*[ r«p*-iPP*lo 

I NO®*  I Of)/ 1 0 

IPO*IPO-INOP*10 

I A la  I NOP Flo 

iNOPalNOP-I AX*10 

IF  IIAI.NE.O)  l A A * l 

IF  IIPO.NE.OT  !PO*I 

IF  (NTn.lT.  2.  OP.  NPM.LT.  2)  I AVP**> 

IF  TIFAO.EO.I*  IAVO*0 

ThFTS*T»*p1 

Ph  I S*  Top? 

OTosT  «P3 
0OmsT«(O4 
Of L0*TNP5 
GNOPaTMPfe 

IFtO»«10 

GO  TO  (4*2.  *.7 .54 .7?  .00)  * I GO 

c 

C EAfCUTF  CAPO  - CALC.  INCLUDING  PAO I ATcO  FIELDS 

C 

38  IF  T IFiOx.FO.IO.ANO.lTHPI .FO.O)  GO  TO  IS 
If  CnFoQ.EO.I.ANO.ITHPI  ,£O.O.ANf'.IFLO*.GT.7T  GO  T l IS 
IF  ( ITvPi .nE .0)  GO  TO  40 
IF  ( I F | Oil.GT.  7|  GO  TO  T9 
!FlO**7 
GO  TO  <.1 
JFLCH*! I 
GO  TO  *1 
IF  AP*0 
PFLOaO. 


*A  ?h\ 

** A ?S? 
Mi  ?*  3 

h A P8<» 

**  A PBS 
PA  ?S© 
Mi  P»7 
MA  PPM 
•*  A p*»9 
MA  ?t  f) 
Mi  291 
PA  ?%? 
Mi  PQ  J 
Mi  PS*. 
*A  PSS 

Mi  PSS 

Mi  ?S7 
HA  ?S8 
HA  ?S9 
HA  TOO 
MA  301 
HA  30? 
HA  103 
MA  10*. 
MA  -»0S 
MA  106 
MA  -*0  7 
MA  108 
MA  109 
MA  11«» 
MA  311 
MA  312 
MA  313 
HA  1J4 
MA  US 
MA  US 
MA  117 
MA  310 
M4  119 
MA  320 
MA  321 
HA  322 
MA  323 

MA  32* 

MA  32S 
MA  32* 
MA  327 
MA  328 
M A 329 
MA  330 
MA  331 
MA  33? 
MA  333 
MA  334 
MA  33S 
MA  336 
MA  337 
MA  338 
MA  339 
MA  340 
MA  341 
MA  342 
MA  343 
MA  344 
MA  345 
MA  346 
MA  347 
MA  34  0 
MA  349 
MA  350 
MA  351 
MA  352 
HA  353 
MA  354 
MA  355 
MA  356 
MA  357 
MA  358 
MA  359 
MA  360 
MA  361 
MA  362 
MA  36  3 
MA  364 
MA  365 
MA  166 
MA  167 
MA  368 
MA  169 
MA  370 
MA  171 
MA  372 
MA  373 
MA  374 
MA  375 
MA  376 
MA  377 
MA  378 
MA  179 


*1 

c 

c 


r 

r 


c 


c 


]>'») 

I A r-*0 

: 3 

1**0 


:t-»» .* 

Iff  »r’-31.*3.?»  P'MSa'.O. 

:»  . * *s : .*£.i»  go  to  *1 

n>-*2 
dp-**9o . 


Mi  ?«2 
Mi  3P3 
Mi  3B* 
Mi  385 
Mi  186 
Mi  ?87 
*•4  i*» 
Mi  189 

Mi  190 

*4  391 
Mi  192 
Mi  193 


GO  ’0  **?. *7. 5<». 72,801  . IGO 
f%3  3*  T«€  “ih  lN«»uT  SECTION 

5Es:v»’*g  re  t*£  f^eq^enc?  do  i oop 


Mi 

M4 

Mi 
Mi 
M 4 


tf  «—7  .E  ..  1 1 GO  TO  MS 

:*■  t Is  = :.e:.i » go  to  mm 
f — z-:Et.f  pq 


Mi 

Mi 

Mi 

Mi 


r.-js-'  — IK^rso 

F5»<-.?/f-w;s 
■l  l»»}*0./r “i? 
p5>'  ;*»S.  Fw^7,4L4M 

d;]h' 

f = Eic€vC*  -CiUsG  OF  GEOMETRIC  oiPiMETEWS 


Mi 

Mi 


Mi 

M4 


- ? ■ s i • *•  Mi 

l ; at  t(»fS  Mi 

T I I > *•  5 ) *ff p Mi 

HU*?'!i»»  Mi 

ST  1 1***1 < :»»F»  M4 

*’•:>*-!  n > *r  p *•* 

!g:*2  Mi 

SEG-ENT  LOinlNG  Mi 

>s;h’  “* 

Ir  C4{.L  f.OiD  (LOTv^.LOTiG.LDTiGF  .L0TiCl*/L9.2Ll.?LC.M  Mi 

i : a : > Mi 

Jff  . M«LGir,E0,»*»  PPJNT  197  Mi 

GPT'.NO  Si  = 4»ffTffP  Mi 

>S:h*  1*9  Mi 

IF  :<S»*p.E0.1»  GO  TO  50  Mi 

Iff  « r=-E=ff.£*.l»  GO  TO  ‘•9  Mi 

yst-Ia-SSP*'  a./(EPS»-S!G*MLAM*S9.Q?*Fj) ) Mi 

:r  ^ilt.-C.O)  GO  TO  *8  Mi 

S:***.*5C-*t.T/»L4,«  Mi 

StS*-»*CP*3T/>L4M  Mi 

®o:hT  j-»9.  NaiDL.SCPwLTtSCffMPT  Mi 

?3>»  1*9  Mi 

53>r  i55.  EOSP.SI6  Mi 

<*-:  •:  =:  Mi 

P= ; * * ;«i  Mi 

r,-  -5  c;  Mi 

P3>T  1*2  Mi 

:>?Ni  *•* 

Mi 

S?-*-CT^E  MiTsji  SET  «»p  Mi 

Mi 


IF  .*  lSTi»T.sC.0»  GO  TO  52  Mi 

IC 1*9  *i 

I C2*l  Cl  *i 

:c3*:ci  m* 

*»t**S.*C*«  Mi 

vCX*VX«  Mi 

Mi 

i Mi 

Mi 

:i...  ITIMli  Mi 

:j«. . :**{ r MPrrf.NCo^cM.NLoiD.PKH) 

C4<._  SF'“»>»0  ITIM2)  Mi 

T;**Ti»2-TI»l  Mi 

Mi 

Mi'PJi  »*CT0»I47TinN  Mi 

Mi 


19* 

195 

396 

*97 

198 

199 
4.00 
*0  1 
*02 
*03 
*0* 
*05 
*06 
*07 
*09 

*09 
*10 
*1  1 
*1? 
*13 
*1* 
*15 
*16 
*17 
*19 
*19 
*20 
*21 
*22 
*23 
*2* 
*25 
*26 
*27 
*28 
*29 
*30 
*31 
*32 
*33 
*3* 
*35 
*36 
*37 
*39 
*39 
**0 
**1 
**? 
**3 
4** 

**5 
**6 
**7 
**0 
**  9 
*50 
*51 

*53 

*5m 

*55 

*66 

*67 


:ic^  f i:TiS<SPEO»V:P*CM*IP.  IX*N30#«NCOL*NrOLS.  IPS’'*'! 


IET1*T*1  Mi  *60 

Jff  ::»•€. LC. II  GO  TO  53  Mi  *M 

s«  ■ « *S-(  i Mi  *6? 

sc;*. *n::c5  ►*  ‘*3 

SEtCNO  ITIMli  Mi  *6* 

T>?*TT«1-TIM?  Mi  *65 

M>T  TJM.TJM2  Mi  *66 

IV* J Mi  *67 

Mi  *69 
Mi  *69 

€*CITar:j%  SET  Up  IPIghT  hand  SIDE.  -E  INC.)  Mi  *70 

Mi  *71 
Mi  *72 

*^wj:*2  Hi  *73 

:sc«;  “a  *?* 

N®*lN**0  Mi  *75 

»r  »:»••?. F0,0»  G*  TO  57  Mi  *76 

Iff  *’P-FL G. lF.O.O*.1*TYP. €0.*)  OPINT  IS*  “i  *T7 

••A«*T|4f»aS  Mi  *70 


" 


copy  r 

FERMI  . . 


HOC  DOES  HOT 

....  „..2E  PRODICTIOH 


TMP4»T  A**PR4 

IF  (IXTYP.lC.I)  00  TO  56 
Tnp|no9|/|tLA<« 

TMP?ai0a?/*LA* 

TSPJaitfQJ/toL*11 

TMPft*«RR6/ ( bLAM*VL A*> 

PRINT  1 Sh . «OR| • 10H?.KPW3. APWu.lPPS, 

00  TO  *7 
fMP|*T««IPR| 

Tmp?»Ta»*pr? 

TMP3*TA*XPR3 

TNPMIPH 

IF  U°7FL0#LC  .*>»  PPINY  155.  YPOl  .jt^.lPJl.hPXI  JlTYP)  ,XP«6 
CALL  Et-NS  iTmpj • T mP7.|nP3vTm04..  I «st  !«| »k.  1<*M • V5ANT . NSANT • 

1CUQ) 

MAf R | K SOLVING  (NfT««  CALLS  S0L*FS1 

IF  (NET.CQ.O.OR.INC.GT.1)  OO  TO  61 

PRINT  ISA 

ITMP3*0 

ITMPliNTYPt l I 

00  GO  |=1.? 

IF  (17*01. FO.  3)  IT-P»=2 
IF  (IfMPl.FO.?)  PRINT  159 
IF  <IT*P1.F0.1»  PRINT  1 GO 

00  59  J«l.NfT 
1Tmp?=nTYP( J) 

IF  ( CIT-P?/ITnp1).FQ.1»  GO  10  S® 

ITMP3af  r*P? 

GO  TO  59 
|TMP4*isEr»l  ( J) 

ITMPS=ISEG?T J» 

IF  (IT»P2.GE.?.AN0.V11 I (jl.LC.G.)  Mil < J>’*LA-*S0OTl<*(!TMPS)-*(IT 
1MP4)  ) ( Y ( |TmpS»-Y  ( T TMP4)  I ••».  ( 1 i | T-aC^*  .7  { ITWP4)  ) **2) 

PRINT  157.  I TAG  ( I TNP4)  • I TmPA, | TAG ( I T--5I  . I f-»5.  YJ  IP  < j>  • T 1 1 ( < Jt  • 

1 V1?R(  J)  .YJ  ?I  C J)  .Y?2R(J)  *Y?2I  ( J)  .»NFT  C2MT-P2-1I  .>NET  (2MTmP2) 

CONT  lNI>E 

IF  ( lT“t>3.FO.O)  GO  TO  tol 

1 TPP1 = I T MPJ 

CONTINUE 

CONTINUE 

IF  ( INC .GT. 1 *ANO. IPTFlG.GT .0)  nopInTbI 

CALL  NET**  CI5FG1.ISEo2.Y11R*yii:.y12p.y1?T.y??o.y??1 »nET  *nTyp* ISA 

Int.vsant. nsant. c*.ip*'*ur.nro*.ncol.ia  .oIH.OlOSNT.nPPInT.NASYM*2P£0 
2.NTS0L . JPMAII 
NT. SOL  = ) 

IF  (IPFO.EQ.O)  GO  TO  62 
ITMP1=mhZ*4»p(mh7-1  > 

IF  (ITmp1.GT.IN09MF-3>>  GO  TO  gp 

FNORHT  TTMP1»»REAL(7PE0| 

FN0RM(ITMPl*l»aAlMAGI7PE01 
FNO RM( fTNP|*?)aCA9SIZPE0) 

FNORM ( t TNP| *3) *C*NG (ZPEO) 

IF  CIPEO.EO.2)  GO  TO  G2 

IF  IFnOPNTITMP1.2>#GT.ZPNORN>  7onCB«*»FvCRm|ITmP1»?) 

CONTINUE 

PRINTING  STRUCTURE  CURRENTS 

IF  f IPTFLG.E0.C-In  GO  TO  64 
IF  < IPTFLG.GT.O)  GO  To  63 
PRINT  161 
PRINT  162 
GO  TO  G« 

IF  (IPTFLG.E0.3.0R.INC.GT.I > GO  TO  64 
PRINT  J63.  APR3.HP0L  CMTYP)  .*006 
PLOSS*(*. 

1 tmp  1*0 
JUMP*IPTFLG*1 
00  70  |*1. H 
CURI *CU®  C 1 1 %#L  AN 
C"AG*CA*S(CU®I > 

PH*C  ANG  f C'JPl  I 

IF  (NL0AD.E0.3)  GO  TO  65 

IF  I ARS  f PEM.  1 7 APR AY  (I))).LT.1.F-?M  GO  TO  G5 
PL05S=PL05S*.5«CNAG»CmAG*PEAL ( 7*R* at (1)1*51(1) 

IF  (JUMP)  69*70.66 
IF  (IPTAG.FO.O)  GO  TO  G7 
IF  (ITAG(I).NF.1PTAG)  GO  TO  >0 
ITMPI*1TMP1*I 

IF  (tTMPl.Lr.IPTAGF.OP. ITMPl.GT, IP7AGT)  GO  TO  70 
IF  ( IPTFLG.EO.O)  GO  TO  69 

IF  ( l PTFLG.LT .2.09. INC »GT .nORmF)  GO  TO  60 

FN0PM(INC)*CNAG 

ISAVEaT 

IF  (IOTFLG.6E.3l  PRINT  164,  «ooi . «PR2 .CmaG.Ph. I 
GO  TO  70 

PRINT  1 65,  I,ITAO(I)  ,*(I)  ,Y(I  ) ,7(1)  .SfU)  ,CUOI  .C**Af,,PH 
CONTINUE 

IF  (MTYP.nE.O)  GO  TO  71 
tmpi*pin-plQnnt-ploss 
TMP2*  1 oO • *r*P!  /PIN 

PRINT  1 G6.  PIN.TMPJ  .pL0S5,PL0GNT,T**p> 

CONTINUE 

IG0*4 

IF  ( I Ft  OV.MF  • 7)  GO  TO  72 

IF  ( l*TYP.6».0. ANO.IXTVP.lt. 4)  GO  TC  11* 

IF  (NFO0.NF.1)  GO  TO  |2l 


39 


s*  ns 

7?  »6 


calculation  of  a • 0 » c in  co*/*«£nt  expansion 


CF.1.  r»aC 
JG*«5 


— ci  co 

rt^  o 
rx?  ~c 

sse  -< 


•v?i»  F rf  1.3  CALCULATION 
IF  (Vfi3.F0.f-l 1 1 00  TO  79 

p»:sr  %*s7 

'ns'*7n*--7np 

TO  T»  1»1.nPZ 

?n«T«7hOT»3ZN0 

:»  ’N£a0.E-).O)  GO  TO  74 

ta*zn-t» 

ST-*5In(T a*ZN0T) 

▼ SOT»f\0*?>fN» 

DC  M J»1.N0V 
»\5T*rS9T*orN0 
IF  IvEaP.CO.O)  GO  TO  TS 
CP-*C0«  tTA*TN»T) 

2B-*S!‘.  <T  4*YN0T) 

«\5’»«sS-*‘»N0 

CO  M «m*1,NPX 
f sCTa«*,^T»0«NO 
JF  'vf»p.EO.O>  GO  TO  76 
i~3*»s=T*5Tm*'‘pm 
r ; s* is; T*5TH»S0H 
ZC®«*nsT»CT« 

GO  TO  77 
■03*XNCT 
rOa*TSPT 
ZOa-Znpt 

T-P|  K*  S/i*l  AM 
*K9>i  ^’3/(1.  AM 

7»»3*7-a/aLAM 

C*Lw  srfL*>  CT-P1  *TmP?.  TNP3,ex»Fv*F.Zl 

T*&I *Ca®S(FX) 

TM5.?*CiHG<FO 

T-P3aCX*S(FY) 

TtoS.aOAsGlEY) 

T*«p;*CJ»SfF2) 

**M»iaCi*»GlFZ| 

FPJNT  168*  *08  • YOB  «Z0B  . T *P1  * T MP2  * T mp3 • T n->w  . T m»S  • " •-•6 

COnTJnvC 

IF  7.EQ.NFP0)  NFAP  = -1 
I f <nF33.nf.i,  g0  t0  79 

8PIN7  135 
*30  TO  15 
COSTIVE 

ST»sCA=0  FaP  FIELD  CALCULATION 

IF  <tF40.EG.-l)  GO  TO  114 
I F UFap.lt. 2)  GO  TO  02 

O0JST 

IF  {IF45.LE.3)  GO  TO  0 1 

: 73,  NOADL.SCB-LT.SC»rf»T 

:r  • :rip.E 3.4)  go  to  «? 

Ir  ’ Is  A3.E3.2.0P.IFAP.E0.S)  MCllFrHPOLdl 
I r < !r •P.E3.3.08. IFAP.EQ.6)  HCtlFsHCIP 
CL*CLT/«LAm 
C***C-T/«LA- 

?3»t:2*C530T  <1./(FPS»2-SIG2**»LAM*5P.0?«F  J)  > 

PP!sT  l?i.  HCLIF.CLI.CHT*EP5P?.SIG2 
IF  <!?is>.NF.l)  GO  TO  ®3 
PPJN7  J75 
GO  TO  aS 
1 

j»i*: 

ITmp* *2*IA»*1 

t’M07«t TMP | .4 

»®IN7  »72  ' 

T*  -F,  o.LT.l .E-20)  GO  TO  04 
F Umi.i,  iH/PFLO 
c I-Aasr l3/bL am 
£I0*2-i*O.ME*»A-MNT  <E*PA»  ) 

*®I*.T  !73.  PFL0.EX88.EXPA 

ppInt  174.  !GTb(II«IGTP(J)«IGax(I  T*»p|  ).1gax(Ivm»3) 

TF  • :«-rP.EG.  )»  GO  TO  07 
IF  »:«*YP.E1.4)  GO  TO  86 
sst:o. 

C*CGN»-.4PI/(  1 .♦*P06»XPP6I 

jZr9*jr CN 
-.-o  it 

P !S«33-.51»X006*X°06*MLAM»41.A8 
GCOP*«*.  AM»«lA»4»2.«dI  / I 376.7  3*PfN» 
p-x: »arN-PLOS5-PLO0NT 
GCON«-,COP 

I r <:°*.NE.O)  GCON»GCON»PIN/PPAO 

I *? 

.cio 

3!n7*7. 

TmP1«03***TA 
t*02*.s*OT«»T4 
o^I»p»«;  s-QPm 
*>0  154  0Pi«l*NPN 


r — 3P» 

-■  ; ra 


40 


pH  I jPh|  »Ow-« 

»T4 

w r*r~f ts-of* 
ro  109  «t*»i*nth 
t *OT«* 

|f  tKS<“P. eO.^.ANO.THrr.iiT. |f4W.NF. I I Gl  t u 109 

T**A*T**«"T*TA 

if  nr»«.ti.n  go  ro  -w 

CALL  FFLO  nH*.5ri*,tT-,|  *-• 

GO  TO  *0 

CALL  GFLO  tOFLO/iL  AM.Oil.T-f  T/ ..  t«(f PH,  ?Q«T  f .<SYNU| 
COOMaCAASCFPOl 
ronA*Ci4Gl£POl 
f THM^r3F4L<ET*4»rONjGIPT«»  » 

F Thh.^PJ  (ctm*?) 

EThAiCanOIFTni 
E Ph1?*u g AL »EP-«CONjG<PMf»  > 

C0HM3SOPT CEPM*2l 

£P-*»CaHG<FPm1 

IF  C I^AP.CO.ll  GO  TO  10* 

ELLIPTICAL  POL API/Al I^S  CALC. 

IF  CE TX-2.GT. 1 ,E-?0.0».  £•>  — ?.  ’•*.!. F-?ni  GO  TO 

TILTA*'). 

r*»AJP?rO. 

FMlS»?rO. 

A«PAT*1. 

I5TNS«i.ALK 

GO  TO  96 

OFA7»E»*A-fTHA 

IF  (FPmA.lT.0.1  GO  TO  92 

0FA22r0FA7-360. 

GO  TO  -31 
0FA22s-.FA2.360. 

IF  C ABSIDFAZ)  .GT.ARSClFAZ?)  i :FiZsr.AA72 
COFA7arOSfOFA7»TA» 

TSTOPi-»eTMM£-et>HM? 

TST0P?=2.*FPHM*FThm»C0FA? 

T IL  T A a . 6*4  TGM?  C TST092  • TST 1P1 I 
STILTAsSlNlTlLTA) 

tsto»i=tsto»i*stilta»st!lta 

TSTO»?9TSTOP?^STlLTA*roSlT!LT»> 

EMAj32,-TST09l*TST0P2*£T-«? 

Em|N«2*TST0O|-TST0O2*pP»’*»? 

|F  CEMfN02.LT .0.1  pmIn92*0. 

A*PAT=^0»T (EMINP2/FMA j32l 

tilta«tilta»to 

I F <AKJA  T.GT.l.E-5)  GO  TO  9- 
ISFNS=*POLC 1) 

GO  TO  96 

IF  (OFA7.GT.O.)  GO  TO  95 
1 SENS  3m»0L 1 21 
GO  TO  96 
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MA 

930 

MA 

931 

MA 

932 

MA 

933 

MA 

936 

MA 

935 

MA 

936 

MA 

937 

MA 

938 

MA 

939 

MA 

960 

MA 

961 

MA 

962 

NA 

96  3 

MA 

966 

MA 

965 

NA 

Q66 

MA 

96  7 

MA 

96  8 

MA 

969 

MA 

950 

MA 

951 

MA 

952 

MA 

953 

MA 

956 

MA 

955 

MA 

956 

MA 

957 

MA 

95* 

MA 

959 

MA 

960 

MA 

061 

MA 

962 

MA 

963 

MA 

966 

MA 

965 

MA 

966 

MA 

967 

MA 

968 

MA 

969 

MA 

970 

MA 

971 

MA 

972 

NA 

973 
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VjafcOiiTTSE  CaOC  CB  1 
C CB  2 
C C*-»C  CA,#P'JTES  COEFFICIENTS  OF  Tm£  CONSTANT  (A).  Sfj£  IB).  AND  CB  ) 
C CSMSt  «C>  T£3HS  IN  T~£  CUWPCNT  iNtf  UPOCAT  {ON  FUNCTIONS  FOB  Th£  f»  * 
C CLIENT  VfCTOO  CO».  CB  S 
C CB  6 


CO^ON  'PAT  A/  N*NP«  A OOO)  • V Mod)  *7  IBOOI  • * l t *•  DO  ) • •»  I I °00  > • *LP  < **<*0  ) . 

1 3fT  » A.**3  » • ] COsl  14301  • |C9N>(800>  . IT  *0  I BO  3 1 #•!.***•  ! »“■  » « • J*ai 
CC^-On  /C«NT/  Al*<SOf»).AI  I I *00 » .B  I * I BOO  * • IMOOl  .C!  «*<*00l  *C  1 1 <*0 
l0>  • Co* I BOO  > 

CO**Os  /Jwnk/  NCO**JOAI2S»  •nCO/«JO/<>S»  .Nf  12*JIZI?S»  CB  u 

CO«?LE*  C'jB.CLO*CLL«CLT  •Ai.flA.r*  C«  12 

DATA  TP/6.2B316S30*/  CB  li 

DO  ?1  I * 1 *N 

call  T3TOU»jC01  •JC02.DIL»DI*1 
ci*"»*“  :l 

CKsTPa- Ik 

S!nl«MS<Cl» 

COSL^SlCl* 

Sts- *S!N«C*I 

co***:  - s «co 

IF  I jcni .OT.I-BOOOO) 10O  TO  s 
jH«*n-( JCOl*BOOOO» 

CLO*C«i*  < J ! * ■<  ) 

GC  TO  ® 

5 IF  ( JCOH  1.6.7 

\ CL9M0..0.) 

|F(s:i f.LT.tlftO  TO  3 

no  2 <*l.sci* 

Jl IKS J* Ilf) 

? CLO»CLr*Cs;B{  JI  AKI 

3 IF  INCO»  ,LT.l»f.O  TO  B 

DO  * *«1«*C0* 

JO  « * » JO  < I *>'  1 

V CL0«CL0-CL'»1  JOA<» 

GO  TO  » 

6 CLO* I 9 • »0. I 
GO  TO  * 

7 CLO«Co»I JCOl » 

IF  « ICON? I JCOl ) .NE . 1 • AnO • jCnl »NF . 1 1CL0*-CLC 
0 CLL  *C’J*  € I A 

IF  A JC02.GT.I-O030OUGO  TO  13 
jIvkss-I JC02*90000) 

cl  » *C‘jp  f ji  a«o 
60  TO  IB 

13  IF  I JC9219.1B.1S 

Q CL**tC..O.) 

IF  INCD.’.LT.IIGO  TO  II 
CO  10  « * l »nC07 
JOZ**JC?<*l 

10  CL^CL^CLPI JOZKI 

11  IF  NCI7.LT. IKO  TO  IB 

CO  12  ■•WnCKZ 

J!?**JT7  CM 

12  CLtsCLT-CuPUIZBl 
GO  TO  IB 

|B  CL**C9.*0.I 

GO  TO  IB 

15  CLt«C-#»«  JC02» 

IFflCONl  (JC021  .NE.I.AN0.jCO2.NF,nCLV»-CLV 

16  IF< JCOl. LT. 1-90000) )GO  TO  17 
IF  I JC32.LT. 1-90000) >GO  TO  1« 

SlL**MNL*COSMCOSL*Sl*K 

CELLO«StNfSlN«-S!L* 

AAs<CLP*5!M*-CLL#SIL«*CLT*BlNL) /CFLLO 

3 As  tCL O*  ICOSB-l . I *CLL« ICOSL-COSW) *CLV*  C I .-COSL  l I /rFLLO 
GO  TO  *3 

17  IF f JC02.LT. 1-90000)  )<>P  TO  19 
SfL*  *COSL*C05k -sinl  *si nk 

cello*mh-cosl 

A«a iCl0*SIN«*CLL*S1L*-CLT»C05l» /CELLO 
B«»lCLO*lCOS«-l.l  *<CLL-ClT»*5lNU /CELLO 
GO  TO  79 

IB  SllxiC^SL^eOSK-SlNL^STN* 

CELL9*r9S*-SlL< 

AisfCLC*COSB-CLL#SlL**CLT#Bl^H.l  /CELLO 

F*  Aa  | l CLL-CLOI  •SINMCLtMI.  -COSL  ) 1 /CELLO 
GO  TO  ?9 

19  CELLO*?. •SISL*fOSL 

AA*CLL* tCL*-CLO>*COSL/CELLO 
BA*ICLT*CL0»»?1BL/CELI.0 

?0  C*aCLL-AA 

A|PIIIa»EAL(AA» 

A11C1»«AINAG1A*I 

0|»(t>*aEAi(A«) 

B1II1»*AI«TAGIBA» 

Cl*MI»««f»LlCA| 

21  ClIlllaAlNAGtCA) 

PE  TjBN 
CNO 
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V7 


SOa&O'JT  !NC  0*5? T tn»0-.NCOt  .C**.**  ' .p«**I  C*  ! 

CM  2 

C**5£T  SETS  T-*E  00*01  E*  STRi^yt-wF  *AJ9J*  jn  T»-E  »=»**  C*  CM  3 

C*«  4 

Common  /G»’ *'  s.sc.x.&001*ri‘i',“>*:,:®j3>.5T(833).affflv*'».  iL°c*00» . 

l 8ftib*oi  .:cc**i  <3*0*  %tcon2<*oo» . me- moo»  .*la*.  ipc  — .j-h 

Di*£*s;0w  C-<**»0»*,»iC0i  » CM  7 

CCmmON  IC»S£  .S?Lr»,,S•s^»t,,  .St  »ST  SrM  CM  8 

com-om  /pcstai/  i:i*::?.ici.*.5f s. •.»-?>•  i-t r-  cm  9 

COMMON  /ANGL'  SAl»(530» 

COMMON  /JUS</  NCO«.^Oxl2S>.A,Cl».jMfPS>.VCOZ.  J07  l ■*«•>  •*•<:!/ .JlZl2S>  CM  11 

COMMON  /P?n./  »H04.PM0Y.P*i0Z.CA3J.9A3j.$ALPO»P*.Pr.:«E?S.S£FP'i  CM  1 ? 

COMMON  //LOAD/  ZAPPA* <800» 

COMMON  /CMC/  7PAl!.ZPATI2*(*L*C-.SCP*I..SCM-tt.NPArL.^Sr-P.IcAM.|Pf»F  C“  14 
0[M£N5fOM  CA9C  | » . ,^A9U  | CM  IS 

OlMENSfOM  £TP13>  . £T  I i 3>  .£?£<?»  .£•*?  «?»  CM  14 

COMPtE*  F i.CM,ZA^oAir.Cl.FPf£T.'Z.*i.C£l.rF?.CE3  CM  17 

COMPLE*  Z3AT!.5tFS*»fFPS.?PS:*..Z3ATiS.Tl,ZSC5N.Z-»’!2  CM  18 

EQUIVALENCE  <CA«1 1 J #A».P<1  ) > • t<A8(!)*4M  1)  »>  .lF2**2r)  . *?»•£«>£>  C«  19 

DATA  ETA/378. ?1/. »!?/**. 283MS3/.P|/?.1-1SP?SS/  CHIP. I 

FUNCTIONS  - - 
JSF  0 C )» * J» ' j-1 ) /np*jP-ax 

JjfQC J) »-j-93000*( <-J- 90001 i/Jt'At.l j *\P 


NOP=M/»*P 
JPMAX*.|M4«/S0P 
J EON*9 

IF  (ICftSE.S*.?)  »E*IV>  11  CM 

FJ*(0..1.»  CM 

!2*2»NPRL***»PO* 

ISLCM*  f 2 Cm 

!T*NP9L*  cm 

MI  1 *ICl *1  CM 

IF  < I C 1 .£C.  0 » 00  TO  I CM 

IF  IICl.LE.C-in  00  T*  2Z  CM 

CALL  aL  C*  I N (1M>I2«1CMI  CM 

CONT  I N<j€  cm 

DO  21  !XSL*1*NII»NPLC^S  C« 

1SV*(  !*9L*1  * l»*NP0L*  CM 

IF  (IiaLKl ,EO.NaLO«S»  I T *Nt  AST  CM 

IF  I IFASE  .LT.  3 ) I T *NCOl  CM 

IF  (N9xCL.E3.«»  CO  10  2 CM 

T1*F  J»?367.067/FLOAT(*.PAOL1  CM 

T2»SCP«P*FlOAT (N®AOL>  CM 

2PATIS*7PATI  CM 

00  3 I * I • I * 

00  3 J*l*sPCM 

CM( J.I)=(0..0,»  CM 

itl*it 

IF  { ISV.CE ,\P) GO  TO  60 

IF  f ISV*IT.OT.SP>ITl*Np-ISV 

CM 

SOU°C€  SC:-m£NT  LOOP  CM 

CM 

00  IB  .1*1  • N CM 

CALL  TOIOl J. JC01* JCC2.0IL*f)IK) 

S*SKJ>  CM 

8*MK  J) 

X J* A ( J)  CM 

VJ*V!J!  CM 

ZJ*2<J>  CM 

CAM J*CAM ( J)  CM 

SAMJ*SA?fJI  CM 

SALPJ*S*L»t J»  CM 

CM 

obsep/ation  segment  LOOP  CM 

CM 

DO  18  TO®  * 1.  ITL  CM 

I* I SV* I po  CM 

nu»*n>**J  cm 

riJMTin-tj  cm 

IJ.I-J  cm 

CAPI*CARIII  CM 

SAp I*S A9  f I » CM 

salpi*salp*!»  cm 

PFL*-1.  CM 

CM 

LOOP  TO  iMTLUr?  IMAGE  OF  SOd®CF  SEGMENT  Fo»  STP*  CTU-t  0*E*  G»OUM>.  CM 

CM 

00  18  !P*1.MS*MP  CM 

PFL*-pfl  CM 

SALPPa^ALPj'Pn.  CM 

Z1  J*ZU»-PrL*FJ  CM 

ZP*XI  )*CA8j*TIJ*SAaj*7| J*SALPP  CM 

01 j*CAP|PCA9j-SAai«SAaj*SAtPI*StLP®  CM 

PMO<*Xf j-CA?J»?P  CM 

PHOT *YI j-SApj*ZP  CM 

PM0Z*ZIJ-SAL»»*7P  CM 

RH*SOPT  |P-C**P“0«»®M0'r#PH0r*P'-O7*P-OZI  CM 

IF  (PN. or. !.€-*»  GO  TO  A CM 

R*»OK*0.  CM 

PHOT*0.  CM 

PmOZ*0.  CM 

Ofp*0.  CM 

GO  TO  a CM 
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21 

23 

24 

25 

26 
2? 
29 

29 

30 

31 
3? 

33 

34 

35 

36 
38 


42 

43 
4* 
*5 

47 

49 

50 

51 

52 

53 

54 

ss 

56 

57 
SB 

59 

60 
61 
62 
S3 
6w 

65 

66 
6? 
68 

69 

70 

71 
>2 

73 

74 

75 

76 

77 

78 

79 

80 
81 
62 
83 
64 
MS 
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ni 
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J** 
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r-1  ijpg' 

— -p* 
CP 
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CP 
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c3 
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1 — 


c 

c 

c 


8 


9 


C 

C 

c 

c 


24 


26 


27 

11 

\2 


?n 

2* 


13 


14 


30 


50 


Uh0Aa9*>0«/P** 

CM 

84 

9* 0»*Jn0</9. 

CM 

87 

9~073®“C//s- 

CM 

PM 

0 I »C  *9  I .9-0<»S4a  I .BHOZ^SAt  PI 

CM 

*9 

P*  SOP  T < ■ I J**  f J*T  I J»  * I J*7IJ*/U) 

CM 

90 

IF  « IP. •.€.2*00  TO  10 

CMRO.l 

I J*1 

CM 

91 

IF  <IpF9F.£?.l»  00  10  10 

CM 

92 

BMAG.B 

CM 

93 

4YM4G*S0PT(Al J**IJ*Tlj»tl J) 

CM 

9<> 

CM 

45 

SET  PARAMETERS  FOP  RAOUL  • !*£  GROUND  SCREEN. 

CM 

56 

CM 

57 

IF  (NP»0l.F3.9>  GO  TO  T 

CM 

93 

XSP£C*»*  « !»•?  J*n  M*Aj>/|H  I»*7JI 

CM 

99 

T5PEC»CTII)*2J*?I!UT jl/l2(1 1 • / II 

CM 

100 

PH0SPC*S0PT«ASPEC*«SPFC*fSPeC*TSP€C*T2*T2> 

CM 

101 

IF  (P-OSPC.0T.SC«m*.  » GO  TO  6 

CM 

102 

2SCRNsTi*B^0SPC*ALCG«PH0SPr/T2l 

CM 

103 

?PATI*» 75CP>«^RAT is)/ (ETA*7PAT IS*/SCPN) 

CM 

134 

GO  TO  7 

CM 

105 

7QAT  W?9»T IS 

CM 

104 

IF  (XYM4G.GT.WE-6!  GO  TO  A 

CM 

107 

CM 

108 

CALCULATION  of  REFLECTION  COEFFICIENTS  «HEN  GpCONC 

IS 

SPECIFIED. 

CM 

109 

CM 

110 

px*o. 

CM 

111 

PT  = 0 • 

CM 

112 

CTm«1. 

CM 

113 

?RSIN»  ( 1 • ••>•> 

CM 

114 

GO  TO  9 

CM 

115 

PX*-Y] j/Xt-40 

CM 

114 

p V *X I J/AVM4G 

CM 

111 

CT-««ZI  J/5M1G 

CM 

113 

2»SIS*:S^»: » 1 ,-ZPAT  U7PAT I*( 1 ,-CTh*CTh) ) 

CM 

119 

REFS«-(CT*—  2»iTI*7»SIn)/(CTh*2R4TI*2»SIN» 

CM 

120 

oft PS* f 7»4TI*CTm-zbsIM /I2»ATI*CTh*2»SIM 

CM 

121 

REFPS»PFro<-j*rs 

CM 

122 

IF (R.GTMa<-»Gn  TO  58 

CM  1 

22.1 

CALL  IST&?A.S.»m./p.OIJ.OIp.ETP.ETI.DIL.OIK.I j. IP. JCO 1 

• jC02.xm.vm, 

7M1 

CM 

124 

fill  M4JRI*  CLEMENTS.  ELEMENT  LOCATIONS  0ET6RMJNF0 

Hr 

CONNECTION 

CM 

125 

OATA. 

CM 

124 

CM 

127 

CCl»C-ot*f£r®ll»*ETI(|») 

CE2*CmPLI (E  TR(2>  *?!l  I?)  I 
CE3*Cmpl*(£T»<3).ETI  (•»>» 

IF  I JCO!  .GT.aijf  |«JSE0(  JCOll 

IF  ( jrf)\  .LT.f-900001  1 J* 1 *JJFQ < JCOl ) 

JE2* JS* Itjl 

IF  f JC02.GT.JI j£3* JSEQ » JC02) 

IF  < JC-V.LT.t-aa0)0) 1 jf3*JJE0<  JC.W 

IF  f JCO!  •'IT.  1-900001  ICO  TO  24 

IFlJC02.LT.  (-90300)100  TO  26 

CE?*CE?-T«/«**CCl 

CE3»CET-I-/I«»-CE1 

CM » JE 1 • I » *Cm  I J£  1 • IPo ) *C£  I /X* 

GO  TO  91 

IF(jCO?.GT.!-QOOOO»»GO  TO  27 

CE2«CE?-t»/»**CF3 

CE1*CE|-7*/im*CCJ 

C»*  I jE3.!pA**CM(  JC3.IP01-CE  3/Xm 

GO  TO  ?7 

0EN*1./  (X»*««4-7M»7M) 

CE2*CE2*CEt-CE3 

CM(Jf  3.I»»»*Cm(JE3.IP»)  ♦0En*(Zm#CE1-XM#CET) 

CMfjfl  .!*»»» »€••(  JE1.IPO)  •0£n*Wm*CE1-?m*CE3) 

IF  f JC01.LT.I-90000HGO  TO  ?9 
IF  € JCCI 111*29.12 

CALL  j-CLS*CEl  .NCIl.Jl*tNC04.  JOX^TPP^CM.NPO^NCOL.NP*  JP**AX| 

GO  TO  '*« 

IFtlCOM?! JtOlI.EO.jlGO  TO  2» 

IF  « JC31.E3.j1G0  To  28 

Cm  I jEI. I »•»*€*•(  JEW  IPR1-CE1 

GO  TO  ?9 

CM i JEI ,1p»*«Cm( jEl.IPo) *CE | 

Cmj  j»2.  f »»»«CM  I JE  ?.i;- I -it  ? 

IF  CJC02.LT.  f-9ftOOO))GO  TO  1M 
!F<jC0?»13.I*.14 

CALL  J-ElSICE  3.NC07*J07*NC1Z* JT7* Ipp»CM*NOO«.NC0L.Np« JPMAXI 
GO  TO  lA 

IFflCO^Ii >C02>  .EO.JlGn  TO  30 
IFf  JC02.C0.J1G0  TO  30 
CM<  if  J.!P®i«CMi  j€3.IPo)-CE3 
GO  TO  1 A 

CM(jCl,li»P»*CM|  JE3.1PR)  *C*  3 


GO  TO  | A C*l44.1 

PA«1«P1?*»  CAU4.2 

A0«7»/*  CM144.3 

Al«»0»f  UWI.-AOPAOH  CAU4.4 

Cl*CM»t  I fCf*5|PnHl  | .-SlMIPAHll ) C**l  44.5 

fO«S*f TA»*f*Cl*C-PL*C  l •«-1./RSm1 I / <P  I 2*P#P»  CHI  44.6 

ET*5»ETA*Al*C|*CMOt Ml  I ..RKHl-l./PKHl I /(2.«p 1 2*P*P)  CM  144. 7 

F7*EP*>0-ET»A)  CMI44.5 

CHI44.9 

|FJ  !P.rQ.2iCAiL  GncE/E  <1WFZ£(?WEp£(1WEPE  (21)  CMJ44.6 

if  ?»  JV9I  J» 

CM|jE2.|>«»»CmoE2,IPo).oI  j*E2*OIP*Ep 


47 


/ jf  >-j 

jrts 

yz  i*  ■ i • s 

j:t l * ! : :*«i  i j» 

I"s2  • .1 

:»  ,C01«900  33.**€  . J€D»»)GO  TO  42 

!*  .C*-?.LT  . <-*D0?9»  ) JCX^mU 
' - -- 

*2  Ir  ' X92*90930.9£.  J£3M»G0  TO  6t» 
X*5»l 

:*  jcoi.lt. «--ooooi » jcas*i 

63  T'-;0<  J.X01  • jC02.DIL.0IK) 

>£H«SlHf2!L»*MSl01«O-SlNt0lL*f>M) 

••»COS’"-lL*C!K) 

! r OI«.lT.OIL)OIL*OIK 

z-«c:s«o:l> 

/:*v 
1 1 /TEN 

?•*  /:e*. 

GO  TO  f47.t*.49.70).JCAS 
67  !*  JC31.ST.3)  Jfi*  JSEQ<  JC01  > 

.**>*  JS^Zi  J» 

^€?. :s>»> *c«t jf?. ipp> *vm/*m 

:•  I *IP®I  •!•/**» 

:r  'XO  1)6-*. 66. 73 


j-fteSfC£!  .NC!<»J!*«NCO**  J0*«  !»9*CM,s.s,C0L  «sp. 


73  :*  re:*.?!  jcoii.e3.jigo  to  71 

!*  jCOI.EO.  J>00  TO  71 

r-  ■ j£i.ra!»)*c-<jEi.fp'>>-2»*/*M 

&:  to  *6 

7 i :•  jEl.!=»9l*C-tJEl.IP£>l*?M/X9 

G3  TO  a* 

*3  -r!  £»-TSEO  < J) 

!*  * jC02.3T.5I J£3*JSE0( JC02> 
r»  ' j€2.  IP®)  *C*I  j£2.!«*0)-ym/<m 
C*  ! .iso|«o»ll.!po»*l./** 

:r  XD?#7%.66.75 
7*  CE3»-7-/*** 

C»«.l  J**ELSCC£?,nC02*  J02*NC!Z*J!7*  IO9*Cm*N.»iC0L*SP. 

TO  4* 

75  !*  •.COsl  f JC02».EO.J>GO  to  7? 

!*  X3P.E0.JI00  To  72 

v€3.:»P)*C-( JE3.IPO) .2M/XH 

SO  TO  65 

72  c*' jCS.TaPJiC-UEJ.IPPJ-ZM/XM 

SC  TO  66 

**  j£?»JS*;tJ) 

.!  ■ X 3 1 • 4 0 0 0 3 1 

0-? ■»*  1 ./  '/“•Z*- **•»»•) 

:•  „£2.i®®»*c«ijE2.sp»)-i. 

C».  I . ::s)arm  T , 1POi-X«4»DEN 

:•  jci.r^oj-c-c  j£i.iP*»)-ZM#o£N 

>0  to  *6 

7J  S*3fJ» 

X :**9-UCC2*40000> 

0£n*  1 . ✓ (Z***7“» **••*••> 

C»  j€Z.!PP)*C«*f  J£?.IPP)*1. 

0»? 1.1=0 )«CM(I.I°P) -K**0£ * 

C«*  ^£!.I®»»*C*»IJEI*1PP)-ZH*0EN 
66  0*T|S>.f 

5*  1*  <»«P*0L.n€.0»  Z»»TI*ZPATIS 

: *i*3I»  ?l£»£*«'5  **OOIFrEO  Pv  LOA^INO 

Ir  'Ml  0*0. EO. 3)  CO  TO  20 

x ;*  i*i. it 

J.ST.SOI&O  TO  23 
:•  :•  j.!)-zabpa<«  j» 

2:  :r  '::a3E.it.3»  go  to  21 

a<. c * 0 • 1 1 1 • 1 1 J2*  1 * 31 ) 

I3:*2*3lM 

:r  '!C1.?9.63L0"S»  GO  TO  ?1 
C-<PPT 

p . ::>*  i**.? 

MCi5E.LT.il  00  TO  22 
9?  • : HO  1 1 
:::«-2 

*r  t JCi^E.EO.S)  IC1»-I 
:-<»oT 

2: 


CM  146 
CM  147 
CM  148 
CM  1 49 
CM  ISO 
CM  1S1 
CM  152 

CM  153 
C"  1*4 
CM  155 
CM  156 
CM  157 
CM  158 
CM  159 
CM  160 
CM  141 
CM  162 
CM  143 
CM  16* 
CM  165 
CM  166- 


«nno 


C 

C 

C 

2 


C 

C 

c 

3 


a 

5 
C 

c 

c 

6 


7 


C 

C 

c 

p 


111 


21? 

?\ 


SU«R0JT1NF  0 A T A >N 

DAT  AGs  IS  T"*t  MAlS  30  jT  IN»  FO®  TNPII?  Of  GEO-ET^r  'ATI. 

COMMON  /©A  | A/  A M3Q>  . * MrtO)  .7<*00»  • ' 1 (BOO)  • »I  ( i'  01  . Ai  PHAOI  • 

1 4*  f l^c  » . | CON | ( 9001  • !C0\2  MOQt  » I T AG  ( 907  I • <L  A***  1 P ■ ’ m.  j».*x 
OlMENSr^s  *2<;  i.  r ^ < 1 1 • 7^(11 
0ImE*»519N  ATST(4>*  |F«<7i»  If  Y ( 2 1 • iF/(2> 
iMTtOC-  gm.atst 

EQUIVALENCE  <<2<1> .St <1)7.  ( i?  I 1 > . *L“  <1 ) ) . tZ2tl»  .«f  T ( n » 

DATA  aTST/2-^#.2mC.a  #7-0W.?mG5,?MGF  . 2hGm • ?» J*  • 2*  jf  / 

DATA  JFA/JW  «/!**  .1HY/.1F7/1H  .1*7/ 

DATA  TA/0.01 *<.532*25?/ 

JPSYN.G 

N*IRE** 

N»0 

PRINT  |1 
POINT  1? 

OEAD  GECMFTOr  DATA  CA^D  anO  SRA.jCh  to  SECTION  FO®  coFRATTQn 
REQUESTED 

READ  1*.  Gm.ITG»sS.**1.y»1.Z«1.xM2.tm2.2»2.RAO 

IF  (GM.FO.AISTU  )>  GO  TO  2 

IF  If.M.to.ATST  (?#  » GC.  TO  3 

IF  (GM.gO. ATST ( 3» I GO  70  - 

IF  (GM.gQ.AT5T  (<.))  GO  TO  5 

IF  (GM.EQ.AT5T (51 ) GO  TO  9 

IF  (Gm.EQ.aTST  (Ml  CO  TO  9 

GO  TO  10 

generate  segment  OAT  A FQR  straight  »I®g. 

NMI«€*Vm!RE*1 
1 1 *N*  1 
I2*N.NS 

PRINT  1-*  N#l®£  » AWl • Y *1 • Z«| • X#?. YW2.7M2*® A0.N5. 1 1 • 12 • I TG 
CALL  MIRE  (A*l .r«1.2»i*A#?#tw2.7M2.0ADtNS.ITG) 

GO  TO  1 

REFLECT  STRUCTURE  ALONG  0®  Z AXES  0®  DOT  ATE  TO  FORM  CvlInOER. 

IY*NS/10 
I Z*NS- I T • I 0 
l X*  I Y ✓ 1 0 
I Y*IY-M*10 
IF  Ua.NE.O)  1**1 
IF  (IV.NE.O)  I Y * 1 
IF  I1Z.nE.0I  17*1 

PRINT  15*  IFX(I**ll.lFY(lY.ll,irZ(IZ*ll*lTG 
GO  TO  - 

PRINT  IS.  NS.ITG 
I4*-l 

CALL  ppflc  (IX.IY,IZ«ITG«NS) 

GO  TO  1 

SCALE  CTRUCTJSE  DIMENSIONS  3y  FACTOR  XMl. 

DO  7 1*1. N 
XtI>*X !l»*X»l 
Y ( I 1 *Y ( I ) *XW 1 

zm*ziii*(«i 

X2lI)a*2(I>M<«l 
V2(Il*r2U)*X-l 
Z2(II*72(I»*»mI 
8I(I)**»I(I)*X-1 
PRINT  1 T . XMl 
GO  TO  ! 

MOVE  STRUCTURE  0®  REPRODUCE  ORIGINAL  STRUCTURE  IN  NE*  POSITIONS. 

PRINT  14.  1TG.sS.spUym1.7m1.XW2.YTi2.Zm2.oaD 

X«1«Xm!*TA 

YM1*YhI»TA 

ZMl«Z«|*TA 

CALL  M«vf  (XAl.YWl.ZMl.*P?.YM2.7W2.1sTtRAO..S) .NS.ITG) 

GO  TO  I 

TERMINATE  STRuCTJOf  GfOm£T»t  In°UT. 

CALL  C^NECT  (ITG» 

CALL  rCNVOT 
jmax*o 

IF (NS.FQ.O)RE  *J»N 
print  on 

REAf)  41  .GM.IX.IV.I7 

print  Q1.GM«IX«1Y,|Z 
IF  (GM.FQ.ATSf  » 7)  10.1  T*>  21 
If  (G*A.s«  .ATST(4)  )GO  Tm  212 
|F  ( jma  » . FQ  # Q | r,Q  TO  20  * 

Of  TU°S 
PRINT  m2 
STOP 

IE (fx.rO.O.ANO.I f.rQ.oiGO  TO  2fl« 
f •ISCGsO 1 1 ■ * I y» 


OA 

1 

OA 

2 

0* 

3 

OA 

CA 

7 

DA 

W 

OA 

9 

DA 

10 

DA 

n 

CA 

12 

OA 

13 

DA 

u 

OA 

is 

CA 

16 

CA 

IT 

OA 

19 

OA 

J9 

OA 

20 

OA 

21 

DA 

22 

CA 

23 

CA 

2A 

DA 

2S 

OA 

26 

OA 

27 

0* 

29 

OA 

2 » 

CA 

30 

OA 

31 

OA 

32 

OA 

33 

OA 

3A 

OA 

35 

CA 

36 

OA 

37 

OA 

39 

OA 

39 

OA 

AO 

CA 

A 1 

OA 

A2 

DA 

A3 

DA 

AA 

OA 

AS 

OA 

A6 

OA 

A? 

CA 

AS 

OA 

A 9 

CA 

SO 

OA 

51 

OA 

52 

OA 

S3 

DA 

Sm 

OA 

55 

DA 

56 

OA 

57 

OA 

50 

OA 

59 

OA 

60 

OA 

61 

OA 

62 

DA 

63 

DA 

6 A 

DA 

65 

OA 

66 

OA 

67 

DA 

69 

OA 

69 

OA 

70 

OA 

71 

DA 

72 

OA 

73 

DA 

7- 

OA 

75 

DA 

76 

OA 

77 

OA 

73 

OA 

79 

OA 

40 

CA 

« l 

DA 

92 

49 


£ 


1 


' 


lr  * I2.‘Q.2)GO  TO  201 
I*  * ICOM  ( I t 

IF  »J«.l  ?.  1-900001  )f.O  TO  1)1 
IF  t l*.fO.O.O»».I*.£O.IlGO  TO  ill 
jMAXsJUAX*) 

12*-  I J*A<  «<»000U  ) 

ICONl ( l I *1 7 
00  TO  204 
201  I » = 1 COM ( I > 

If (I*. LT. <-90000)100  "0  111 

rf ( ix.Fo.o.otf.M.fo.iioo  »o  in 

J"H*  J«AX*1 

r/*-<j-*x.o3oooi 
IC0N21 1 1*17 

204  !r  « !«.»  T.0100  TO  2*5 

if  » icom  c iu  .ea.iiGJ  to  ic*> 
lCCN2im*I  2 
00  TO  III 

206  ! CONI  1 1 « l * I Z 
GO  To  111 

205  00  207  1*1. N 

If  < ICOM  < I)  .tQ.m  ICOnI  (I  >*IZ 
If  r IC0**2<  ft. CO.  (OlCOs2(  1**12 

207  COST  ISi’C 
GO  TD  HI 

208  If  ( J*A«  .61.0)00  TO  212 
•00  210  1*1.  N 

IX*0 
r Y=-I 

t2s-tJ*4X*9009l) 

00  20Q  j* 1 »N 

!f  (ICOsl  (JI.NF.IYIGO  TO  211 
1**1 

TCCNl ( Jl * 1 7 
GO  TO  209 

211  If  MC0N2UI  .NE.IT1G0  200 
1**1 

ICOM  < j) *1 Z 
200  CONTINUE 

If  Cl*. fO. 0)00  TO  213 
jmax*j*ax»1 
2io  continue 

213  If (GM.fQ.ATST (81 1PFTU9N 
GO  TO  1)1 


10  PRINT  10  0*  84 

P3 1 NT  20«  r»M.  I TG»N«..X*1  .Y*l  .2«1  .Xw2.Vto2.Z*?.OAO  0*  65 

STOP  0*  06 

C DA  87 

|l  FO»MAT  f////*.33X»3SH STPJCTliaf  SPEC  If  IC*T  ION  - - -.//.37x.  OA  88 

1 20HCCO»OINATES  MUST  8£  INPUT  IN./.37*.  DA  89 

2 29MMETEPS  OP  BE  SCAlEO  TO  NET£PS. / . 37X*  OA  90 

3 31M9EF0PE  STRUCTURE  iNPvJT  IS  FNOEC./Fl  OA  91 

CODMAT  ( 2*.«»HMIPE*79*.8MN0.  Of  . AX  ,SNF  I&ST  *2*  • A-'.  • ST  . 5 * • 3*T  AG*  0*  92 

1 /•2*.THN0..a*.2M<l.‘>*«2Hri,9*.2H71»liM.2HX2*9*.2-T2.v«.2-22»6X.  DA  93 

2 S*OA3IuS.3X.‘»HSEG.»5».4MSEG..‘»X*<*»iSFG*.5<.lHN0.>  DA  94 

13  FOOMAT  (A2. 13.ISWf  10.5)  DA  95 

14  fOaMAT  (lx.I5«3fll.5*lX,4f 11.5.9X.T5.4X.I5.|X«15«*>X*15)  OA  96 

15  fOPMAT  C6*»34HSTauCTUPf  PEfLECTEO  ALONG  THF  AXES*"1 1 l X • A l ) • DA  97 

1 2?h.  TAGS  INC^EmentEo  By. IS)  DA  98 

16  fCPMAT  (6X . 30mSI4uCTu-E  ROTATED  ABOUT  Z-AXIS.I3.  DA  99 

1 JO*  TT-ES.  LA8LES  KCpEmFnTED  BY, IS)  DA  100 

17  FORMAT  C6X.26MSTauCTuat  SCalED  BY  F ACTOR • f 1 0 .S ) OA  101 

18  FOPMAT  (6X.49HTHE  STRUCTURE  HAS  PEf  N movEO.  MOVE  ~ATA  ^APO  IS  -/  DA  102 

l 6**I3. I5.7F10.5)  OA  103 

19  FOPMAT  (25m  GEOMET»Y  DATA  CAPO  FPPOP)  DA  104 

20  FOPMAT  (lX.A2.I3«I5«7rlo.5)  OA  105 

90  FOPMAT (////. 46H  JUNCTIONS  USING  OPT  tONAL.  INT EPPCL * T I On  *ETH)0> 

91  fOPMAT (A2.I3.2lS) 

93  FOPMAT ( lx. A2, 13.215) 


92  FOPMAT (69H  £RPOP—  INVALID  OATA  CAPO  MNEHE  JUNCTION  INTERPOLATION 
1CAP0  F»P£CTEO) 

ENO  OA  106- 


50 


ei 

u 

El 

II 

D 

n 

D 

[I 

El 

El 

[I 

Ei 

D 

n 

B 

b 

D 


i 


copy  avail;.:’/  t”  pp5  eoes  not 

PERMIT  FULLY  L'T^l  . LGJUCTION 


mm 


Su*4  O'JUnE  F*CT«  <s.  A. 

SU-*ODJTlsF  TO  FACT  - t •*'-]*  T*jT*>  A *JN  I T To  J ANGi  iL  Arf  «ATPl* 

ASO  as  'jPPfB  T»I  As'-^c  A-  S!*.r-  T~E  GAiiSS-30«L!fTLF  ALOO-IT^i 

p=fsemfd  on  page s *1  :-*!*•  ^ a.  3*lsTon--a  fiw,t  course  in 

Stj«C»KAl  ANA,  tsis.  -ClO-  ^FFE*  to  COMxrtgTS  jN  RALSTONS 

If  * 7 • <*»A7wr«  T-iss-Jj;'  . 

0 1 MENTION  A I NO I •< «N?  1 **  » • NCt-4, 

FA 
F A 
FA 
F A 
FA 
FA 
FA 
FA 
FA 

\ 

2 

3 

4 

5 

6 
7 
A 
5 

CO«***ON  ZSCoATm/ 

CO-PLf*  A.O.APJ 

f A 

1! 

IsTEC.Fo  P.ohi  .PP1  j.ps 

FA 

12 

jr,  f.,o 

FA 

13 

OQ  0 B=l.N 

FA 

1- 

FA 

15 

ST?p  1 

FA 

16 

F A 

17 

1 <*|.N 

F A 

:* 

C")=1'P.K| 

F A 

i? 

COnTIV.F 

F A 

23 

FA 

21 

STfPS  3 ASO  3 

FA 

22 

F A 

23 

9*.  1 a Q — | 

FA 

?•. 

Ir  <**«KLT.l>  GO  T r - 

FA 

25 

CO  3 

F A 

26 

Pj* 1° ‘ J) 

FA 

27 

AO JrO(Pj| 

FA 

23 

A <P. J) sAQJ 

0(PJ»=0( J1 

FA 

23 

Jal*J*I 

FA 

30 

00  2 I s J® 1 .N 

FA 

31 

0U*»0[  I )-AU.n*t;  J 

cost  im.f 

FA 

33 

CONTINUE 

FA 

34 

COnT  Iv  E 

F A 

35 

FA 

36 

STEP  - 

FA 

37 

FA 

35 

OuAixPfAL  <0<R)  •C0n~j<*>|3)  ) > 

FA 

33 

IP|»»  *s 

FA 

40 

**»!*»•  1 

FA 

41 

IF  IPPKGT.NI  GO  TO 

F A 

42 

DO  5 l*»Pl*S 

FA 

43 

ELMAG*3EALI0U)*CCn^G  «:;!)!» 

FA 

44 

IF  <ELNAG.LT. 0*aAA)  C*C  *2  5 

FA 

45 

OMAUEl.NAG 

FA 

46 

IP  IP) *1 

FA 

47 

CONTINUE 

FA 

48 

COST  IV, E 

FA 

49 

IF  lONAA.LT.l.C-lOl  JFLG»1 

FA 

50 

99s 1 P ( 9 | 

FA 

51 

A (O.P,  :f)(PO) 

FA 

52 

0<99)sO(P> 

FA 

53 

- A 

54 

STFP  s’ 

FA 

55 

FA 

56 

IF  <PPt.GT.N)  GO  TO  • 

FA 

57 

A3J*K/A<P.B» 

DO  7 I»0PJ. N 

FA 

50 

A(p.I)30<I)#A«J 

FA 

55 

continue 

FA 

60 

CONTINUE 

FA 

61 

IF  <IFLG.EO.Ol  GO  TO  * 

FA 

62 

PPfNT  10*  P.OMAA 

FA 

63 

IFLG*3 

FA 

64 

CONTINUE 

FA 

65 

P£Tl  l°N 

• 

FA 

66 

FA 

67 

C/)OM  at  IlH  «6HPIV‘* 

FA 

6* 

CNO 

FA 

69- 

SURROUT  INE  lNrr.(8,S.RM.7P.0I  J.OIR.ETO.ET  l .OIl.OK  , J !»•  jCOl . JC02 

1 


i 


C INTO  COMPUTES  Three  C^hplEs  FlflO  COMPONENT*;  T*.»T  «*u. T»pli  The 

C ThbEE  SEGMENT  CURRENTS  USfQ  IN  J N TF  SPOL  A T I S'*,  0/?B  * Sfj-£NT. 

C The  SE  COMPONENTS,  t T p u NO  E T I • r-0  INTO  Th£  JME-aMJA  MATRIX. 

c 

OImENSTON  FTR<3>.  ETIlll 
DATA  10/6.283135308/ 

C 

C CO-PUTF  TANGENTMt  FIELD  On  03^FhvAT(0N  SEGMENT  O f TO  SINE. 

C COSINE.  ANO  CONSTANT  '*U»R£nTS  On  SO-jOCF  SES^EST. 

C 

CALL  EFLD  IP.S.MM.?p, t J.E/PS.E7TS.Eq»S.EPTS.E7DC.r 2 IC.E59C.EOIC. 
1 E?R<.F?I*c.ER»S.EBt* I 
IF  (IP.NE.2)  GO  TO  l 
CALL  Gn  (E7RS.F?IS.£PoS.E*'IS) 

CALL  Gn  (E7PC.E7IC.EWOC.F-ICI 
Call  Gn  (f  7HA.f7K.£w=S.f  RfAJ 

1 ETos*E7PS*0! J*E»RS*OlR 
£T!S«EMS*DI J.Fols.OIa 
FT»C»E^RC»OI J*ERMC*OIo 
CTIC*E7IC*0I  J.£«IC*0I*» 

E Tok>£7  5k *ni J*ERR«*Olo 
f TJrf*E*I**0!  J*EPI^*0I3 

C 

C COMPUTE  INTERPOLATION  COEFFICIENTS  AnO  F03m  T h£  COEFFICIENTS  OF 

C The  ThpEE  SEGMENT  CURRENTS  uSEO  IN  CURRENT  INTERPOLATION. 

C 

CL*TP*0IL  * 

C«*TP*DI< 

SINL*SIN«CL» 

C0SL*COS<CL» 

SlNKsSfNICMI 

cosK-crs(CM) 

SILK»SIN<CL*C<1 

COnS*SINL*SIN*-SIl»< 

ETO(1)*(SINK*FT»'<*TCOSS-1.)*ETDS-SINK*ETRC)/CONS 
ETMI)*ISINK*ETIS*(C0SK-1.)*ETIS-SINK»ETIC)/C0NS 
ETR(2>*(-S!LK*ETR<.(C0SL-C0SK)*ETRS*  (S|NL*SfN»O»£*iO/C0NS 
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1 
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Integer  r.ri .r2.pr.s-  lf 

LOGICAL  Ll.ia.L3  LF 

IFlC»*0  if 

C LF 

C INITIAI  ire  »I.P?.;l.J>  LF 
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IF  UCAsE.fO.  T.OR.ICAsE.EG.S)  r.f\  7 ->  R 
CALL  SOLVE  (NPEQ.A( Irt.i).lP(lA)  . * ( ! A > . s-0. ) 
00  10  10 

CALL  LFSOLV  ( A •N«,rO*NCOL  .IjUlA'.^t’Al  ) 
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IF  (N3P.fQ.ll  GO  TO 
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00  11  «*l*NOP 
IA*I*f«-l»«SPFQ 
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S*  5 1 T J)  10 
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